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Abstract. Rigorous pointwise asymptotics are established for semiclassical soliton ensembles (SSEs) of 
the focusing nonlinear Schrodinger equation using techniques of asymptotic analysis of matrix Riemann- 
Hilbert problems. The accumulation of poles in the eigenfunction is handled using a new method in which 
the residues are simultaneously interpolated at the poles by two distinct interpolants. The results justify 
the WKB approximation for the nonself adjoint Zakharov-Shabat operator with real-analytic, bell-shaped, 
even potentials. The new technique introduced in this paper is applicable to other problems as well: (i) it 
can be used to provide a unified treatment by Riemann-Hilbert methods of the zero-dispersion limit of the 
Korteweg-de Vries equation with negative (soliton generating) initial data as studied by Lax, Levermore, and 
Venakides, and (ii) it allows one to compute rigorous strong asymptotics for systems of discrete orthogonal 
polynomials. 



1. Introduction 

Many important problems in the theory of integrable systems and approximation theory can be recast as 
Riemann-Hilbert problems for a matrix- valued unknown. Via the connection with approximation theory, and 
specifically the theory of orthogonal polynomials, one can also study problems from the theory of random 
matrix ensembles and combinatorics. Roughly speaking, solving a Riemann-Hilbert problem amounts to 
reconstructing a sectionally meromorphic matrix from given homogeneous multiplicative "jump conditions" 
at the boundary contours of the domains of meromorphy, from "principal part data" given at the prescribed 
singularities, and from a normalization condition. So, many asymptotic questions in integrable systems {e.g. 
long time behavior and singular perturbation theory) and approximation theory {e.g. behavior of orthogonal 
polynomials in the limit of large degree) amount to determining asymptotic properties of the solution matrix 
of a Riemann-Hilbert problem from given asymptotics of the jump conditions and principal part data. 

In recent years a collection of techniques has emerged for studying certain asymptotic problems of this 
sort. These techniques are analogous to familiar asymptotic methods for expanding oscillatory integrals, 
and we often refer to them as "steepest-descent" methods. The basic method first appeared in the work of 



Deift and Zhou |DZ93|. The first applications were to Riemann-Hilbert problems without poles, in which 
the solution matrix is sectionally holomorphic. Later, some problems were studied in which there were a 
number of poles — a number held fixed in the limit of interest — in the solution matrix (see, for example, the 



paper |DKKZ96| on the long-time behavior of the Toda lattice with rarefaction initial data). The previous 
methods were extended to these more complicated problems through the device of making a local change 
of variable near each pole in some small domain containing the pole. The change of variable is chosen so 
that it has the effect of removing the pole at the cost of introducing an explicit jump on the boundary of 
the domain around the pole in which the transformation is made. The result is a Riemann-Hilbert problem 
for a sectionally holomorphic matrix, which can be solved asymptotically by pre-existing "steepest-descent" 
methods. Recovery of an approximation for the original sectionally meromorphic matrix unknown involves 
putting back the poles by reversing the explicit change of variables that was designed to get rid of them to 
begin with. 

Yet another category of Riemann-Hilbert problems consists of those problems where the number of poles 
is not fixed, but becomes large in the limit of interest, with the poles accumulating on some closed set 



F in the finite complex plane. A problem of this sort has been addressed | KMMOC | by making an explicit 
transformation of the type described above in a single fixed domain G that contains the locus of accumulation 
F of all the poles. The transformation is chosen to get rid of all the poles at once. In order to specify it. 
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discrete data related to the residues of the poles must be interpolated at the corresponding poles by a 
function that is analytic and nonvanishing in all of G. Once the poles have been removed in this way, the 
Riemann-Hilbert problem becomes one for a sectionally holomorphic matrix, with a jump at the boundary 
of G given in terms of the explicit change of variables. In this way, the poles are "swept out" from F to the 
boundary of G resulting in an analytic jump. There is a strong analogy in this procedure with the concept 
of balayage (meaning "sweeping out" ) from potential theory. 

In establishing asymptotic formulae for such Riemann-Hilbert problems, it is essential that one make 
judicious use of the freedom to place the boundary of the domain in which one removes the poles from 
the problem. Placing this boundary contour in the correct position in the complex plane allows one to 
convert oscillations into exponential decay in such a way that the errors in the asymptotics can be rigorously 
controlled. If the poles accumulate with some smooth density on F C G, the characterization of the correct 
location of the boundary of G can be determined by first passing to a continuum limit of the pole distribution 
in the resulting jump matrix on the boundary of G, and then applying analytic techniques or variational 
methods. The continuum limit is justified as long as the boundary of G remains separated from F. 

This idea leads to an interesting question. What happens if the boundary of G as determined from passing 
to the continuum limit turns out to intersect F7 Far from being a hypothetical possibility, this situation is 
known to occur in at least three different problems: 

1. The semiclassica l limit o f the focusing nonlinear Schrodinger hierarchy with decaying 
initial data. See ICMMOOt . This is an inverse-scattering problem for the nonself adjoint Zakharov- 



Shabat operator. On an ad hoc basis, one replaces the true spectral data for the given initial condition 
with a formal WKB approximation. There is no jump in the Riemann-Hilbert problem associated with 
inverse-scattering for the modified spectral data, but there are poles accumulating asymptotically with 
the WKB density of states on an interval F of the imaginary axis in the complex plane of the eigenvalue. 
The methods described above turn out to yield rigorous asymptotics for this modified inverse-scattering 
problem as long as the independent time variable in the equation is not zero. For t = 0, the argument 
of passing to the continuum limit in the pole density leads one to choose the boundary of G to coincide 
in part with the interval F. Strangely, if one sets t = in the problem from the beginning, an 



alternative method due to Lax and Levermore |LL83] and extended to the nonselfadjoint Zakharov 



Shabat operator with real potentials by Ercolani, Jin, Levermore, and MacEvoy [EJLM93| can be 
used to carry out the asymptotic analysis in this special case; this alternative method is not based on 
matrix Riemann-Hilbert problems, and therefore when taken together with the methods described in 



|KMMOO' does not result in a uniform treatment of the semiclassical limit for all x and t. At the same 



time, the Lax-Levermore method that applies when t = fails in this problem when t ^ 0. 
The zero-dispersion limit of the Korteweg-de Vries equation with potential well initial 



data. As pointed out above, the original treatment of this problem by Lax and Levermore [LL83] was 
not based on asymptotic analysis for a matrix-valued Riemann-Hilbert problem. But it is possible to 
pose the inverse-scattering problem with modified (WKB) spectral data as a matrix- valued Riemann- 
Hilbert problem and ask whether the "steepest descent" techniques for such problems could be used 
to reproduce and/or strengthen the original asymptotic results of Lax and Levermore. In particular, 
we might point out that the Lax-Levermore method only gives weak limits of the conserved densities, 



and that a modification due to Venakides [V9C| is required to extract any pointwise asymptotics {i.e. 
to reconstruct the microstruture of the modulated and rapidly oscillatory wavetrains giving rise to the 
leading-order weak asymptotics). On the other hand, "steepest descent" techniques for matrix- valued 
Riemann-Hilbert problems typically give pointwise asymptotics automatically. It would therefore be 
most useful if these techniques could be applied to provide a new and unified approach to this problem. 

If one tries to enclose the locus of accumulation of poles (WKB eigenvalues for the Schrodinger 
operator with a potential well) with a contour and determine the optimal location of this contour for 
zero-dispersion asymptotics, it turns out that the contour must contain the support of a certain weighted 
logarithmic equilibrium measure. It is a well-known consequence of the Lax-Levermore theory that the 
support of this measure is a subset of the interval of accumulation of WKB eigenvalues. Consequently, 
the enclosing contour "wants" to lie right on top of the poles in this problem, and the approach fails. 
In a sense this failure of the "steepest descent" method is more serious than in the analogous problem 
for the focusing nonlinear Schrodinger equation because the contour is in the wrong place for all values 
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of X and t (the independent variables of the problem) , whereas in the focusing nonlinear Schrodinger 
problem the method fails generically only for t = 0. 

The large degree limit of certain systems of discrete orthogonal polynomials. Fokas, Its, and 



Kitaev |FIK92| have shown that the problem of reconstructing the orthogonal polynomials associated 
with a given continuous weight function can be expressed as a matrix-valued Riemann-Hilbert problem. 
It is not difficult to modify their construction to the case when the weight function is a sum of Dirac 
masses. The corresponding matrix- valued Riemann-Hilbert problem has no jump, but has poles at 
the support nodes of the weight. The solution of this Riemann-Hilbert problem gives in this case the 
associated family of discrete orthogonal polynomials. If one takes the nodes of support of the discrete 
weight to be distributed asymptotically in some systematic way, then it is natural to ask whether 
"steepest descent" methods applied to the corresponding Riemann-Hilbert problem with poles could 
yield accurate asymptotic formulae for the discrete orthogonal polynomials in the limit of large degree. 



Indeed, similar asymptotics were obtained in the continuous weight case [ DKMVZ99 using precisely 
these methods. 

Unfortunately, when the poles are encircled and the optimal contour is sought, it turns out again 
to be necessary that the contour contain the support of a certain weighted logarithmic equilibrium 



measure (see [ KR98 for a description of this measure) which is supported on a subset of the interval 
of accumulation of the nodes of orthogonalization (i.e., the poles). For this reason, the method based 
on matrix-valued Riemann-Hilbert problems would appear to fail. 

In this paper, we present a new technique in the theory of "steepest descent" asymptotic analysis for 
matrix Riemann-Hilbert problems that solves all three problems mentioned above in a general framework. 
We illustrate the method in detail for the first case described above: the inverse-scattering problem for the 
nonself adjoint Zakharov-Shabat operator with modified (WKB) spectral data, which amounts to a treatment 
of the semiclassical limit for the focusing nonlinear Schrodinger equation at the initial instant t — Q. This 
work thus fills in a gap in the arguments in KMMOO] connecting the rigorous asymptotic analysis carried 



out there with the initial-value problem for the focusing nonlinear Schrodinger equation. Application of the 
same techniques to the zero dispersion limit of the Korteweg-de Vries equation will be the topic of a future 
paper, and a study of asymptotics for discrete orthogonal polynomials using these methods is already in 



preparation PKMMOl 



The initial- value problem for the focusing nonlinear Schrodinger equation is 



subject to the initial condition ip{x,{)) — ipoix). In |KMMOC|, this problem is considered for cases when the 



initial data ipo{x) = A{x) where A{x) is some positive real function R — > (0, A]. The function A{x) is taken 
to decay rapidly at infinity and to be even in x with a single genuine maximum at a; = 0. Thus A(0) = A, 
A'{0) = 0, and A"{0) < 0. Also, the function A{x) is taken to be real-analytic. With this given initial data, 
one has a unique solution of (|l]) for each h > 0. To study the semiclassical limit then means determining 
asymptotic properties of the family of solutions 'il'{x, t) as ?i | 0. 

This problem is associated with the scattering and inverse-scattering theory for the nonselfadjoint Zakharov- 
Shabat eigenvalue problem |^S72 |: 



h— = —iXu 4- Aixjv 
dx 

(2) 

h— = —A(x)u + iXv 
dx 

for auxiliary functions u[x\X) and v{x]X). The complex number A is a spectral parameter. Under the 
conditions on A{x) described above, it is known only that for each H > Q the discrete spectrum of this 
problem is invariant under complex conjugation and refiection through the origin. However, a formal WKB 
method applied to (^ suggests for small h a distribution of eigenvalues that are confined to the imaginary 
axis. The same method suggests that the reflection coefficient for scattering states obtained for real A is 
small beyond all orders. 
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It is therefore natural to propose a modification of the problem. Rather than studying the inverse- 
scattering problem for the true spectral data (which is not known) , simply replace the true spectral data by 
its formal WKB approximation in which the eigenvalues are given by a quantization rule of Bohr-Sommerfeld 
type, and in which the reflection coefficient is neglected entirely. For each h> 0, this modified spectral data 
is the true spectral data for some other (?i-dependent) initial condition V'o(^)- Since there is no reflection 
coefficient in the modified problem, it turns out that for each h the solution of (|l|) corresponding to the 
modified initial data ipo{x) is an exact iV-soliton solution, with N ^ h'~^. We call such a family of A^-soliton 
solutions, all obtained from the same function A{x) by a WKB procedure, a semiclassical soliton ensemble, 



or SSE for short. We will be more precise about this idea in § g. In |KMMOC|, the asymptotic behavior of 
SSEs was studied for t ^ Q. Although the results were rigorous, it was not possible to deduce anything about 
the true initial- value problem for (|l|) with ipo{x) = A{x) because the asymptotic method failed for t = 0. In 
this paper, we will explain the following new result. 

Theorem 1. Let A{x) be real- analytic, even, and decaying with a single genuine maximum at x — 0. Let 
iPq{x) be for each h > the exact initial value of the SSE corresponding to A{x) (see § Then, there exists 
a sequence of values of h, h — for N — 1,2,3,..., such that 

(3) lim Hn ^0 

N^oo 

and such that for all x 0, there exists a constant > such that 

(4) \i;'^-{x)-A{x)\<K.,hf-\ for iV = 1,2,3,... 
for all v > 0. 

As iPq{x) is obtained by an inverse-scattering procedure applied to WKB spectral data, this theorem 
establishes in a sense the validity of the WKB approximation for the Zakharov-Shabat eigenvalue problem 
(^. It says that the true spectral data and the formally approximate spectral data generate, via inverse- 
scattering, potentials in the Zakharov-Shabat problem that are pointwise close. The omission of a; = is 
merely technical; a procedure slightly different from that we will explain in this paper is needed to handle 
this special case. We will indicate as we proceed the modifications that are necessary to extend the result 
to the whole real line. The pointwise nature of the asymptotics is important; the variational methods used 



in [ EJLM93 suggest convergence only in the L sense. Rigorous statements about the nature of the WKB 
approximation for the Zakharov-Shabat problem are especially significant because the operator in is 
nonself adjoint and the spectrum is not confined to any axis; furthermore Sturm-Liouville oscillation theory 
does not apply. 



2. Characterization of SSEs 



Each iV-soliton solution of the focusing nonlinear Schrodinger equation (|I|) can be found as the solution 
of a meromorphic Rieniann-Hilbert problem with no jumps; that is, a problem whose solution matrix is a 
rational function of A G C. The A^-soliton solution depends on a set of discrete data. Given N complex 
numbers Ao,...,AAr_i in the upper half-plane (these turn out to be discrete eigenvalues of the spectral 
problem (^), and N nonzero constants 7o,...,7jv-i (which turn out to be related to auxiliary discrete 
spectrum for (^), and an index J = ±1, one considers the matrix m(A) solving the following problem: 

Riemann-Hilbert Problem 1 (Meromorphic problem). Find a matrix m(A) with the following two prop- 
erties: 

1. Rationality: m(A) is a rational function of X, with simple poles confined to the values {Xk} o-nd the 
complex conjugates. At the singularities: 



Res m(A) = lim m(A)CTi 

A=Afc A^A. 



(5) 



Res m(A) 

A=Ar. 



lim m(A)(Ti 

A^AT ^ 





Ck{x,t) 

-ck{x,ty 
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(6) 



for k — 0, . . . , N ~ 1, with 



N-\ 
^ \ n=0 



Ik 



2. Normalization: 



(7) 



N-l 



■exp{2iJ{XkX + \lt)/h) 



m(A) ^ I, 



Here, cti denotes one of the Pauli matrices: 



(8) 



" 


1 ■ 




■ -i' 




■ 1 " 


1 





0-2 


i 


0-3 := 


-1 



The function 'ip{x, t) defined from m(A) by the Umit 

(9) ■>P{x,t) = 2i lim Atoi2(A) 

A — >oo 

is the iV-sohton solution of the focusing nonlinear Schrodinger equation (|l|) corresponding to the data {A^} 
and {7fe}. 

The index J will be present throughout this work, so it is worth explaining its role from the start. It turns 
out that if J = +1, then the solution m(A) of Riemann-Hilbert Problem |^ has the property that for all fixed 
A distinct from the poles, m(A) — > I as a; ^ +oo. Likewise, if J = — 1, then m(A) ^ I a,s x ^ — oo. So as 
far as scattering theory is concerned, the index J indicates an arbitrary choice of whether we are performing 
scattering "from the right" or "from the left" . Both versions of scattering theory yield the same function 
tp{x,t) via the relation (|9|), and are in this sense equivalent. However, the inverse-scattering problem involves 
the independent variables x and t for (|l|) as parameters, and it may be the case that for different choices of 
X and t, different choices of the parameter J may be more convenient for asymptotic analysis of the matrix 
m(A) solv ing Riemann-Hilbert Problem |l|. That this is indeed the case was observed and documented in 
i KMMOOt . So we need the freedom to choose the index J, and therefore we need to carry it along in our 
calculations. 

A semiclassical soliton ensemble (SSE) is a family of particular iV-soliton solutions of (|^) indexed by 
N — 1,2,3,4,... that are formally associated with given initial data V'o{2;) = A(x) via an ad hoc WKB 
approximation of the spectrum of (^. Note that the initial data ipoix) = A(x) may not exactly correspond 
to a pure A^-soliton solution of (Q) for any h, and similarly that typically none of the iV-soliton solutions 
making up the SSE associated with ipo{x) = A{x) will agree with this given initial data at t = 0. 

We will now describe the discrete data {A^} and {^k} that generate, via the solution of Riemann-Hilbert 
Problem ^ and the subsequent use of the formula (H), the SSE associated with a function ^po{x) = A(x). 
We suppose that A{x) is an even function of x that has a single maximum at a; = 0, and is therefore "bell- 
shaped" . We will need ^(a;) to be rapidly decreasing for large x, and we will suppose that the maximum 
A := A{0) is genuine in that A"{Q) < 0. Most importantly in what follows, we will assume that A{x) is a 
real-analytic function of x. 

The starting point is the definition of the WKB eigenvalue density function p^{rj): 



(10) 



TT 



dx 



-iv) 



y/A{x)'^ + rf 



defined for positive imaginary numbers rj in the interval (0, iA), where X- (77) and x+{'q) are the (unique by our 
assumptions) negative and positive values of x for which iA{x) = rj. The WKB eigenvalues asymptotically 
fill out the interval (0, iA), and ^"(77) is their asymptotic density. This function inherits analyticity properties 
in rj from those of A(x) via the functions x±{'r]). Our assumption that A{x) is real-analytic makes p°(?7) 
an analytic function of rj in its imaginary interval of definition. Also, our assumption that A{x) should be 
rapidly decreasing makes p^ir]) analytic at 77 = 0, and our assumption that A{x) have nonvanishing curvature 
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at its maximum makes p^{r]) analytic at 77 = lA. From this function it is convenient to define a measure of 
the number of WKB eigenvalues between a point A G (0, iA) on the imaginary axis and lA: 

(11) e°{X):=-n r p%rj)dv. 



J\ 

Now, each A^-sohton solution in the SSE for A{x) wiU be associated with a particular value h — hj^, 
namely 

1 1 

(12) n^hN:=--j^ P'^^'^^'^'^^N^J ^(^^"^^ 

where N E Z+ . In this sense we are taking the values of H themselves to be "quantized" . Clearly for any 
given A{x), = 0{1/N) which goes to zero as N becomes large. For each N E Z+, we then define the 
WKB eigenvalues formally associated with A{x) according to the Bohr-Sommerfeld rule 

(13) 0°{Xk) ^T:hN{k + l/2), for fc = 0, 1, 2, . . . , TV - 1 
and the auxiliary scattering data by 

V ,^K ( »(2K+l)g°(Afc) 

(14) 7fc := ~i{-lr exp 

V n-N 

Here, K is an arbitrary integer. Clearly the Bohr-Sommerfeld rule ( |l3|) implies that choosing different integer 
values of K in ( |l^ will yield the same set of numbers {7^}. However, we take the point of view that the 
right-hand side of (p^ ) furnishes an analytic function that interpolates the {7fc} at the {A^-}; for different 
K E 1 these are different interpolating functions which is a freedom that we will exploit to our advantage. 
In fact, we will only need to consider = or A' = — 1. 

For A(x) given as above, the SSE is a sequence of exact solutions of (|l|) such that the A^th element 
'0^"(a;,i) of the SSE (i) solves (|l|) with h = Hn as given by ( [l^ ) and (ii) is defined as the A^-soliton solution 



corresponding to the eigenvalues {Afc} given by (13) and the auxiliary spectrum {7fc} given by ( |l4| ) via the 
solution of Riemann-Hilbert Problem |] with h — Hn. For each N, we restrict the SSE to t = to obtain 
functions 

(15) Vo"(x) :=^''"(a;,0). 

It is this sequence of functions that is the subject of Theorem ^ In the following sections we will set up a 
new framework for the asymptotic analysis of SSEs in the limit N 00, a problem closely related to the 
computation of asymptotics of solutions of (|l|) for fixed initial data ipoi^) = ^i^) in the semiclassical limit. 

3. Removal of the Poles 

The asymptotic method we will now develop for studying Riemann-Hilbert Problem ^ f or SSEs is especially 
well-adapted to studying the case of i = 0, where the method described in detail in p<MMOC | fails. To 



illustrate the new method, we therefore set t = in the rest of this paper. Also, we anticipate the utility of 
tying the value of the parameter J — ±1 to the remaining independent variable x by setting 

(16) J :— sign(a;) . 

In all subsequent formulae in which the index J appears it should be assumed to be assigned a definite value 
according to (|l6|). 

We now want to convert Riemann-Hilbert Problem |^ into a new Riemann-Hilbert problem for a sec- 
tionally holomorphic matrix so that the "steepest-descent" methods can be applied. As mentioned in the 
introduction, in KMMOO| this transformation can be accomplished by encircling the locus of accumulation 



of the poles, here the imaginary interval (0, lA), with a loop contour in the upper half-plane and making a 
specific change of variables based on the interpolation formula (|lj) for some value of if € Z in the interior 
of the region enclosed by the loop and also in the complex-conjugate region. One then tries to choose the 
position of the loop contour in the complex plane that is best adapted to asymptotic analysis of the resulting 
holomorphic Riemann-Hilbert problem. The trouble with this approach is that it turns out that for t = the 
"correct" placement of the contour requires that part of it should lie on a subset of the imaginary interval 
(0, ij4), that is, right on top of the accumulating poles! For such a choice of the loop contour, the boundary 
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values taken by the transformed matrix on the outside of the loop would be singular and the "steepest 
descent" theory would not apply. 

So taking the point of view that making any particular choice of A' G Z in ( p^ ) leads to problems, we 
propose to simultaneously make use of two distinct values of K in passing to a Riemann-Hilbert problem 
for a sectionally holomorphic matrix. Consider the contours illustrated in Figure 0, arranged such that 




Figure 1. The geometry of contours introduced in the complex X-plane. The uppermost 
common point of the contours C'l Cm and Cr is X = iA. The six-fold self-intersection 
point is the origin A = 0. 



{Ao, . . . , Aat-i} C DlU Dr. For A € Dl, set 
(17) M(A) m{X)a^ 




cxp 



2iX\x\-ie"{X)\ 



IN 



For A e Dr, set 



(18) 



M(A) := m(A)CTi 




exp 



2iX\x\+ie°{X) 



IN 



1 



For Xe DlUD*j^ set M(A) := (T2M(A*)*(T2, and for all other complex A set M(A) = m(A). So rather than 
enclosing the poles in a loop and making a single change of variables inside, we are splitting the region inside 
the loop in half, and we arc using dilfcrcnt intcrpolants ([l4| ) of the {7*:} at the {Afe} in each half of the loop. 
Some of the properties of the transformed matrix M(A) arc the following. 

Proposition 1. The matrix ]V[(A) is analytic m C \ S where S is the union of the contours Cl, Cr, and 
Cm, and their complex conjugates. Moreover, M(A) takes continuous boundary values on E. 
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Proof: The function 6'° (A) is analytic in Dl and Dj^ if Cl and Cr are chosen close enough to the 
imaginary axis since ^'^(77) is analytic there. By using the residue relation (^) and the interpolation formula 
( p^ ) alternatively for ii' = and K — one checks directly that the poles of m(A) are canceled by the 
explicit Blaschke factors in and (|8|). □ 



Proposition 2. Let M±(A) denote the boundary values taken on the oriented contour Yi, where the subscript 
"+" (respectively "—") indicates the boundary value taken from the left (respectively from the right). Then 
forXeCL, 



(19) 



For A e Cr, 



M_(A)-iM+(A) = 







n 

\k=0 



X-Xk 



exp 



2zA|a;| - z6>"(A) 



(20) 



For XeCM, 



M_(A)-iM+(A) = CTi^ 




exp 



2iX\x\+i9°{X) 



(21) 



M_(A)-iM+(A) = a^' 




exp 



2jA|x| 



I'M 



■ 2 cos 



0«(A) 



'TV 



On t/ie contours in the lower half-plane the jump relations are determined by the symmetry M(A) — 
CT2M(A*)*cr2. All jump matrices are analytic functions in the vicinity of their respective contours. 

Proof: This is also a direct consequence of (|l^) and (^8|). The analyticity is clear on Cl and Cr since 
0^{X) is analytic there, while on Cm one observes that as a consequence of the Bohr-Sommcrfeld quantization 
condition (|l3|), the cosine factor precisely cancels the poles on Cm contributed by the product of Blaschke 
factors. □ 

Although we have specified the contour Cm to coincide with a segment of the imaginary axis, the reader 
will see that the same statements concerning the analyticity of M(A) and the continuity of the boundary 
values on S also hold when Cm is taken to be absolutely any smooth contour in the upper half-plane 
connecting A = to A = iA. Given a choice of Cm, the contours Cl and Cr must be such that the topology 
of Figure]^ is preserved. We also have specified that Cl and Cr should lie sufficiently close to Cm (a distance 
independent of ?ijv) so that ^"(A) is analytic in Dl and Dr. Later we will also exploit the proximity of these 
two contours to Cm to deduce decay properties of certain analytic functions on these contours from their 
oscillation properties on Cm by the Cauchy-Riemann equations. 

Taken together. Proposition |^ and Proposition |^ indicate that the matrix M(A) satisfies a Riemann- 
Hilbert problem without poles, but instead having explicit homogeneous jump relations on S given by the 
matrix functions on the right-hand sides of (|l9|), ([20|), and (pT]). The normalization of M(A) at infinity is 
the same as that of m(A) since no transformation has been made outside a compact set, so if M(A) can be 
recovered from its jump relations and normalization condition, then the SSE itself can be obtained for t = 
from (||) with m(A) replaced by M(A). 



4. The Complex Phase Function 



We now introduce a further change of dependent variable involving a scalar function that is meant to 
capture the dominant asymptotics for the problem. Let 5(A) be a complex-valued function that is independent 
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of H, analytic for A e C \ (Cm U C|/) taking continuous boundary values, satisfies g{X) + g{X*)* = 0, and 
5(00) = 0. Setting 

(22) N(A) := M(A) cxp(- g{X)<T3/n) 

we find that for A G Cl, 



(23) 

where 

(24) 



N_(A)-iN+(A)-a, 



1 

aL(A) 1 



'1 ) 



ai(A) i 




exp 



/ 2iA|a:| -i6'°(A)-2Jg(A) \ 



Similarly, for A G Cr, we find 
(25) N_(A)-iN+(A) 
where 
(26) 

Finally, for A G Cm, 



1 

aniX) 1 



'1 



V, A;=0 



X-Xu. 



exp 



2iX\x\+ie°{X) - 2Jff(A) \ 



(27) 

where 

(28) 

and 
(29) 



N_(A)-iN+(A) = ai 



exp(i6'(A)/ftAr) 

aM(A) exp(-i6'(A)/?iiv) 



/Af-l 



■ I TJ - K\ f 2iX\x\ - Jg+^X) - Jg-{X)\ ^ fe°{X) 
aniX) * ( 11 x-\l j ^''P ' I • 2cos ' ^ ' 



In 



9{X) :=zJ(.9+(A)-g_(A)). 



This means that given a function g(X) with the properties described above, one finds that the matrix 
N(A) satisfies another holomorphic Riemann-Hilbert problem with jump conditions determined from ([2^), 
(pq), and ( |27| ). Because g(oo) = and g(A) is analytic near infinity, it follows that the correct normalization 
condition for N(A) is again that N(A) I as A ^ 00. These same conditions on g{X) show that if N(A) 
can be found from its jump conditions and normalization condition, then the SSE can be found via (^ with 
m(A) replaced by N(A). 

The function g{X) is called a complex phase function. The advantage of introducing it into the problem 
is that by choosing it correctly, the jump matrices (p^), (25), and ( |27| ) can be cast into a form that is 
especially convenient for analysis in the semiclassical limit of — > 0. The idea of introducing the complex 
phase function to assist in finding the leading-order asymptotics and controlling the error in this way first 



appeared in the paper |DVZ94] as a modification of the "steepest-descent" method proposed in |DZ93| 



5. PoiNTwiSE Semiclassical Asymptotics of the Jump Matrices 

For our purposes, we would like to have each element of the jump matrix for N(A) of the form cxp(/(A) //ijv) 
for some appropriate function /(A) that is independent of ft^v. While this is not true strictly speaking, it 
becomes a good approximation in the limit Hn with A held fixed (the approximation is not uniform 
near A = or A = iA). In this section, we describe the pointwise asymptotics of the jump matrix for N(A) 
with the aim of writing all nonzero matrix elements asymptotically in the form exp(/(A)/?iAr) with a small 
relative error whose magnitude we can estimate. 

Roughly speaking, the intuition is that the product over k of Blaschke factors should be replaced with 
an exponential of a sum over k of logarithms. The latter sum goes over to an integral that scales like hjj^ 
in the semiclassical limit. On the contour Cm, the cosine that cancels the poles must also be encorporated 
into the asymptotics. 
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The branch of the logarithm that is convenient to use here is most conveniently viewed as a function of 
two complex variables: 

(30) L'^,{X)■.^\ogi~^iX-fJ)) + '-^. 

As a function of A for fixed ?/, it is a logarithm that is cut downwards in the negative imaginary direction from 
the logarithmic pole at A = 77. Equivalently, i^(A) can be viewed as the branch of the multivalued function 
log(A — 77) for which arg(A — 77) G (— 7r/2, 37r/2). Suppose 77 G Cm- The boundary value of i^(A) taken on 
Cm as A approaches from the left (respectively right) side is denoted by L^_^{X) (respectively L°_(A)). The 

average of these two boundary values is denoted by i^j(A). 

All the results we need will come from studying the asymptotic behavior of two quotients: 



(31) ^(A) := (jl exp (^-^ (^J^^ L',Wp"{v) dv + f L^.WWY dv 



and 
(32) 



, fe=0 

The function S{X) is analytic and nonvanishing for A G C+ \ Cm- Let us denote hy Q d C+ the domain of 
analyticity of p^{X) restricted to the upper half-plane, so that by our assumptions on A{x), Cm C fl. Then, 
due to the zeros of the cosine on the imaginary axis, which match the poles of the product below X = iA 
and arc not cancelled above A = iA, T(A) is analytic and nonvanishing for A € $7 \ y, where V is the vertical 
ray from A = to infinity along the positive imaginary axis. The domain of analyticity for T(A) is a subset 
of Q rather than of the whole upper half-plane due to the presence of the averages of the logarithms in the 
integrand of (|3^). Whereas these are boundary values defined a priori only on Cm, the integrals extend 
from Cm to analytic functions in the domain \ V via the introduction of the function ^'^(A) (c/. equation 

Lemma 1. For all X in the upper half-plane with hjf < |9?(A)| < B, where B is positive and sufficiently 
small, but fixed as Hn 0, 



(33) S'(A) = 1 + 



[mx)\j 



Proof: Let us define the function 771(7;) by 



777(77) :=- / p'iOdC 



(34) 

This analytic function takes the imaginary interval [0, iA] to the real interval [0, M] where 

(35) M = m{iA) = - / A{x) dx . 

I" J -00 

Since ^"(C) does not vanish on Cm, we have the inverse function 77 = e(777) defined for 771 near the real 
interval [0, M]. Using these tools, we get the following representation for 5(A): 



(36) 

and 
(37) 



S{X) = exp(-/(A)) where /(A) = ^ /fe(A) , 



fc=0 



4(A) := 



1 



mfc + fijv/2 



e(m)(A) - -t'e(m)(A) 



dm — 



-^-e(mfc)(A) — ie(mfc)(A) 



ASYMPTOTICS OF SEMICLASSICAL SOLITON ENSEMBLES: RIGOROUS JUSTIFICATION OF WKB 



11 



with rrik '■= M — h]\[{k + 1/2). Expanding the logarithms, we find that 



1 



(38) Ik{\) = ^ / dm dC I 



2e"(0A3 - 2e"(0e(0'A + 4e'(e)2e(0A' 



(A2 - 



This quantity is clearly 0{h\) for A fixed away from Cm- Now, when |5R(A)| — o(l) as Kn I 0, we can 
estimate the denominator in the integrand to obtain two different bounds: 

(39) 2e'-(0A^ - 2e"(g)e(0^A + 4e'(0^e(e)A _ ^ ( 1 



(A2-e(^)2)2 \^{\f 

and 

2e"(0A3-2e-(e)e(02A + 4e-(02e(e)A ^ / 1 
^ ^ (A2- 6(0^)2 l^|^3(A)-e(^)|^ 

The idea is to use the estimate ( ^9[ ) when e{mk) is close to z3(A) and to use the estimate ( ^0| ) for the 
remaining terms. Suppose first 3(A) is bounded between and A, i.e. there are small fixed positive numbers 
5i and 62 so that 5i < 3(A) < A — S2, and let e = e(/iAr) be a small positive scale tied to h and satisfying 
fiN ^ e ^ 1, and let Li be chosen from 0, . . . , — 1 so that e(rnLj) is as close as possible to i(3(A) + e), 
and likewise let L2 be chosen from 0, . . . , — 1 so that e{mL.^) is as close as possible to i(3(A) — e). Using 
( ^ ) we then find that 

(41) Y.h{\) = 

k=Li 

because the sum contains 0{€/hiq) terms and the volume of the region of integration for each term is 0(?i^), 
and we must take into account the overall factor of I/Hn. Now in each of the remaining terms /fc(A), we 
have 

^^^^ |z3(A)-e(0P " ^ ((mfe-m(z3(A)))^ 

so using ( ^ ) and summing over k we get both 

(43) ^4(A)^0(^^) and ^ /,(A) = O ( ^ 

The total estimate of /(A) is then optimized by a dominant balance among the three partial sums. This 
balance requires taking e ~ |5R(A)|, upon which we deduce that under our assumptions on A, we indeed have 



3fi(A)2 



(44) ^(^) = «^^) and consequently ^(A) - 1 ^ O (^) , 

when 3(A) is bounded between and A. When 3(A) w or 3(A) ~ A, the estimate ( ^9| ) should be used 
only for those terms that correspond to m near zero or m near M respectively. In both of these exceptional 
cases, the same estimate is found. When 3(A) is bounded below by A, there is no need to use the estimate 
( ^ ) at all, and the relative error is of order Sjv uniformly in 3?(A). This completes the proof. □ 

We now use this information about S{X) to effectively replace the sums of logarithms by integrals, at least 
on some portions of the contour S. 

Proposition 3. Suppose that the contour Cl is independent of Hn and that for some sufficiently small 
positive number B, Cl lies in the strip ~B < 3fi(A) < and meets the imaginary axis only at its endpoints 
and does so transversely. Then 

HA rO W 



(45) 



OLiX) = zexp \^l-\^2zX\x\+ L^^iX) {^) dr^ + J ^Ll{X)p^ {r,*)* dr^ - 2J g{X) 



\X\) \\X-iA\ 



as Hn goes to zero through positive values, for all X G Cl with \X\ > Hn and \X — iA\ > Hn. 
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Proposition 4. Suppose that the contour Cr is independent of Hn and that for some sufficiently small 
positive number B, Cr lies in the strip < 5R(A) < B and meets the imaginary axis only at its endpoints 
and does so transversely. Then 



auiX) = *cxp (-^(2*A|.t|+ r^L«(A)pO(7y)d77+ r L^^i\)p'> (rj*)* drj ~ 2Jg{\) 



(46) 




iA 



fiN J \ \\X\J \\X-iA\^ 
as Hn goes to zero through positive values, for all X Cz Cr with \X\ > Hn and \X ~ iA\ > Hn . 

Proof of Propositions [| and ^; These propositions follow directly from Lemma |l] upon using the transver- 
sality of the intersections with the imaginary axis to replace 0(1/|5R(A)|) by 0(1/|A|) + 0(1/|A — iA\). □ 

We notice that the first factor on the second line in ( p5| ) and the first factor on the second line in ( ^ ) 
are both exponentially small as Hn goes to zero through positive values, as a consequence of the fact that 
p'^irj) drj is an analytic negative real measure on Cm ■ This follows from the Cauchy-Riemann equations and 
the geometry of Figure It will be a very useful fact for us shortly. 

Now we turn our attention to the function T(A). The result analogous to Lemma |^ is the following. 

Lemma 2. For all A in the upper half-plane with < |5R(A)| < B, where B is positive and sufficiently 
small, but fixed as Hn — > 0, 



(47) T(A) = l + 



\mx)\ 



Proof: We begin with the jump condition 



(48) 



L%{X)p"{v)drj+ f L%{X)p°{vrdv 

iA 



Ll_{X)p\r^)dr,+ f i°_(A)p"(r?*)*d,7-2»0°(A). 



-iA 

relating the boundary values of the logarithm i"(A) on the imaginary axis. Using this jump relation and the 

definition of i°(A) as the average of the boundary values of L°_,_(A) and L"_(A), we see that for 3?(A) < 0, 
we have 

2ie^(X) 



(49) T(A) = Six) ( 1 + exp 
while for 5R(A) > 0, we have 

(50) T(A) = S{X) ( 1 + exp 



2^g°(A) 

Now, using the fact that ^^(77) is an analytic function satisfying (0^(77) G iM+ for rj E Cm, we see by the 
Cauchy-Riemann equations that in both cases, the exponential relative error term is of the order e"^!'*^'*')!/''™ 
for some K > 0. Since this is negligible compared with the relative error associated with the asymptotic 
approximation of S{X) given in Lemma ^ the proof is complete. □ 

Unfortunately, we need asymptotic information about T(A) right on the imaginary axis, which contains 
the contour Cm , so we need to improve upon Lemma ||. We begin to extract this additional information by 
noting that under some circumstances, it is easy to show that T(A) remains bounded in the vicinity of the 
imaginary axis. 

Lemma 3. // either (i) X is real or (ii) \X\ = A and 3(A) > 0, and if for some B > sufficiently small 
|5R(A)| < B, then T(X) is uniformly bounded as Hn — > 0. 
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Proof: It suffices to show that S{X) is bounded under the same assumptions, because from (^9|) and ( |50| ) 
and the Cauchy-Riemann equations, we see easily that |r(A)| < 2|S'(A)|. 
Using the function m(-) and its inverse e(-), we have the following 

7V-1 

(51) 

fc=0 

where 

A + e(m) 



hN log \S[X)\ = H{mk)hN - / H[m) dm , 

j,-n "'0 



(52) 



H(rn) := log 



A — e(m) 

0, and therefore |S'(A)| = 1 and hence |r(A)| < 2. 



When A e K, we see immediately that H{m) 

Now consider A — iAe^^ with 6 sufficiently small independent of h^f. The idea is that of the terms on the 
right-hand side of (jsj), the discrete sum is a Riemann sum approximation to the integral. The Riemann 
sum is constructed using the midpoints of N equal subintervals as sample points. If H"{m) is bounded 
uniformly, then this sort of Riemann sum provides an approximation to the integral that is of order N""^ or 
equivalently In this case, we deduce that S{\) = 1 + OIHn) and in particular this is bounded as ?ijv 
tends to zero. But as A approaches the imaginary axis, the accuracy of the approximation is lost. 

For A ~ iAe^^ , the function Him) satisfies H{Q) = H'{M) = and takes its maximum when m = M, 
with a maximum value 

(53) 7f(M) = log|cot(6'/2)| . 

Therefore, as tends to zero, H{m) becomes unbounded, growing logarithmically in 0. As a consequence 
of this blowup the approximation of the integral by the Riemann sum based on midpoints for |A| = A fails 
to be second-order accurate uniformly in 0. However, because the maximum of H(rn) always occurs at the 
right endpoint, it is easy to see that when the error becomes larger than 0(?i^) in magnitude its sign is 
such that the Riemann sum is always an underestimate of the value of the integral, and consequently the 
right-hand side of (|l|) is negative. This is concretely illustrated in Figure || where we have taken the example 



N=10, Gaussian Data 



N = 20, Gaussian Data 



(3=0.0125 

0=0.05 

0=0.2 



0=0,0125 

0=0.05 

0=0.2 



Figure 2. The midpoint rule Riemann sums approximating the integral, pictured here for 
the Gaussian initial data A{x) = \pKe~^ . When the peak of H(m) becomes underresolved 
for small 6, the Riemann sums underestimate the value of the integral by an amount that is 
of the order h^. 

of the Gaussian function A{x) — \pKer^ in order to supply the function p^ij]) and therefore the function 
e(m) needed to build H(m). In this case, A — ^/tt and M — 1. The error of the Riemann sum is worst 
when = 0. In this case it is easy to see that the discrepancy contributed by only the subinterval adjacent 
to the logarithmic singularity of H{m) is (1 — log2)fi,Ar + 0(/i^), which clearly dominates the 0{h^) error 
contributed by the majority of the subintervals bounded away from m = M. Consequently, for those A on 
the circle |A| = A for which log |>5'(A)| is not asymptotically small in Hn, it is negative, and therefore S{X) is 
uniformly bounded for |A| = A, as is T(A). □ 
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Using this information, we can finally extract enough information about T(A) on the imaginary axis to 
approximate aA/(A) for A G Cm- 

Proposition 5. Let Cm be a fixed contour from X = Q to \ = lA lying between Cl and Cr, possibly 
coinciding with the imaginary axis. Then, for ^ > arbitrarily small, 



(54) 



1 / riA M 



omW = ic^p\j- [ 2zX\x\ + I L;iX)p%rj) + j <(A)p°(7y*)* dr; - Jg+(A) - Jg^iX) 



■N \ Jo J-iA 



|A| ; ' \\X~iA\)) ' 
as goes to zero through positive values, for all X G Cm with \X\ > and \X — iA\ > hjq . 

Proof: Let C be the closed contour illustrated in Figure ^. This counter-clockwise oriented contour 



lA 



C 



' ' 



Figure 3. The contour C of the Cauchy integral argument. 



consists of two vertical segments, one horizontal segment that lies on the real axis, and an arc of the circle 
of radius A centered at the origin. The function r(A) is analytic on the interior of C and is continuous on 
C itself. In fact it is analytic on most of the boundary, failing to be analytic only at A = and A ~ iA. 
Therefore for any A in the interior, we may write 

If we let Cin denote the part of C with |3?(s)| < Hn, and let Cout denote the remaining portion of C, then 
we get 



(56) 
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Using the estimate guaranteed by Lemma |^ in the integral over Cout, and the uniform boundcdness of T(s) 
(and therefore of T(s) — 1) guaranteed by Lemma ^ in the integral over C-m, we find 

(57) 



\T{X) - l\ < Ki,,hN sup 

seCin \s - M 



KoutflN^OghN sup - 

seCout \s - A| 



for some positive constants Ki^ and if out- Replacing the logarithm by a slightly cruder estimate of ft^^ for 
arbitrarily small positive ^ completes the proof. □ 

We have therefore succeeded in showing that, at least away from the self- intersection points of the contour 
E, the jump matrices for N(A) as defined by ( p3[ ) for \ £ Cl, ( p5| ) for A G Cr, and (|2^) for A G Cm are 
well-approximated in the semiclassical limit Hn ^ hy matrices in which all nonzero matrix elements are of 
the form exp(/(A)/S7v) with /(A) being independent of Hn- The fact that this approximation is valid even 
when the "active" contour Cm is taken to be right on top of t he poles o f the meromorphic Riemann-Hilbert 
problem for m(A) is an advantage over the approach taken in | KMMOC |. 

Using these approximations, we can introduce an ad hoc approximation of the matrix N(A). First, define 



lA 



(58) ^(A) := 2*A|a:| + / 4(A)p"(ry) ^77 + / ^i;(A)p°(77*)* ^77 - Jg+(A) - Jg-(A) , for 



A € C 



M ; 



and for A e Cl or Cr, define 

(59) r(A) := 2a|a;| + 

Then we pose the following problem 



i°(A)p°(77)d7y+ / L0(A)p"(r;*)*d7y-2J5(A). 



-iA 



Riemann-Hilbert Problem 2 (Formal Continuum Limit). Given a complex phase function g{X) find a 
matrix N(A) satisfying: 

1. Analyticity: N(A) is analytic for A G C \ S. 

2. Boundary behavior: N(A) assumes continuous boundary values on E. 

3. Jump conditions: The boundary values taken on E satisfy 



(60) 



N+(A)-N_(A)r7i 







I exp 



T{x)-ie"{x) 



IN 



for Xe Cl, 



(61) 



N+(A)=N_(A)ai 



I exp 



-(A) + 100(A) 



for X G Cr, and 



(62) 



N+(A)=N_(A)ai 



exp 



I exp 



z0(A)A 
h-N J 



exp 



i9{X) 



for X G Cm- For all other A G E (that is, in the lower half-plane), the jump is determined by the 
symmetry N(A) = cr2N(A*)*(T2. 
4. Normalization: N(A) is normalized at infinity: 



(63) 



N(A) ^ I as X 
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6. Choosing g{X) to Arrive at an Outer Model 

Let R{X) be defined by the equation i?(A)^ = + the fact that R{X) is an analytic function for A 

away from the imaginary interval / [—iA{x),iA{x)], and the normalization that for large A, i?(A) ~ —A. 
For ?7 G In Cm, let 

(64) pi,):=Av) + ^'^ ' ^-^^^ ^"(^^'^ 



iv - s)R{s) ni J,a^^) (t] - s)R{s) ' 

It is easy to check directly that for all G / n Cm, we have p(r/) e iR+. Also, using the fact that p'^{s) 
is purely imaginary on the imaginary axis, and that R{s) is purely imaginary in the domain of integration, 
where it satisfies R{—s) = —R{s), we see that 

(65) P(0)=p"(0). 
Furthermore, it follows easily from ( |6^ ) that for all ry G / n Cj\/, we have 

(66) < -ip{rj) < -ip\r,) , 

with the lower constraint being achieved only at the endpoint^ of /, A = iA[x)^ and the upper constraint 
being achieved only at the origin in accordance with (|65|). 
Now, set 

(67) 5(A) := ^ f Ll{\)p{i^*r dr^ + ^ f^^^^ L^iX)p{r,) drj . 

^ J-iA(x) ^ Jo 

This function satisfies all of the basic criteria set out earlier: it is analytic in C \ {Cm U C^.j) and takes 
continuous boundary values, it satisfies g{X) + g{X*)* = 0, and it satisfies g{oo) — because 



Piv*)*dv+ / piv)dTj = 0. 

iA{x) Jo 

Note that ^(A) is analytic across Cm for A above iA{x). Consequently 9{X) = for all such A. For A G Cm 
below iA{x), 9{X) becomes (c/. equation (p9[)) 

(69) 0{X) = -7T / p{v) dTJ . 



We now describe a number of important consequences of our choice of 5(A). 
Proposition 6. For all X e I n Cm = [0, iA{x)], 0(A) = 0. 
To prove the proposition, we first point out that 

(70) lim 4>(X) = , 

AGCa/ 

simply as a consequence of the fact that both p^{r]) and p{'r]) are purely imaginary on Cm- Next we point 
out that 

(71) 0'(A)=O 

whenever A G [0, iA{x)]. This follows from a direct calculation in which all integrals are evaluated by residues 
and the formula ( p^ is used. 

Next we consider 4'{X) for A G Cm \ [0,iA{x)], that is, above the endpoint of the support. Clearly, 
(/)(A) + i9{X) is the boundary value on Cm of an analytic function defined near Cj\/ in D^. Since the 
boundary value taken below the endpoint is i9{X) because 0(A) = there, and the boundary value taken 
above the endpoint is 0(A) because 9{X) = there, we obtain the formula 

piA{x) 

(72) ^{X) ^ i9+{X) ^ ~tn p+{v)dv 



^It is often convenient to think of the function p{ri) being extended to all of Cm by setting p{r]) = for A above the endpoint 
iA(x). In this case one views the lower constraint as being active on the whole imaginary interval lA]. 
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valid for A S Cm above iA{x), where by p+{ri) for 77 in the imaginary interval {iA{x),iA) we mean the function 
p{vi) defined by ( |6^ for rj in the imaginary interval (0, analytically continued from (0,i^(x)) in the 

clockwise direction about the endpoint A — iA{x). In particular for such A we have 

(73) ^'(A)=z^p+(A). 

Carrying out the analytic continuation, we find from that for ry G {iA{x)^ lA), 



(74) 



P+(A) 



R{\) pO(g*)*^g i?(A) 

»A (A - S)i?(s) TTi 



P.V. 



p°(s)ds 
(-^ - s)R{s) 



From this formula we see easily that for all A strictly above the endpoint iA{x)^ p+{X) is positive real. 
Consequently, from ( |73| ) and since 0(A) = for A = iA{x), we get the following result. 

Proposition 7. The function 4>{X) is negative real and decreasing in the positive imaginary direction for 
AeCM\[0,iA(x)]. 

Now we consider the behavior of the function t(A) on Cl and Cr. From the definitions of the functions 
t(A) and 0(A), we see that for A G Cl, 



(75) 

and for A G Cr, 
(76) 



t(A) = 4>(X)+i0(X)~i0°{X), 



t{x) = 4>{x) - ie{x) + ie'^ix) . 



That is, the analytic function t(A) takes boundary values from the left on Cm equal to (j>{X) + i9{X) — i9^{X) 
and from the right on Cm equal to 0(A) — i9{X) + i9^{X). First consider the situation to the left or right of 
the imaginary interval [0,i^(a;)]. Since 0(A) = in [0,ij4(x)], the function t(A) on C^ will be the analytic 
continuation of i9{X) — ?0°(A) from Cm and the function r(A) on Cji will be the analytic continuation of 
—i9{X) + i0°(A) from Cm- From ( [6^ ) we see that for rj € [0, z>l(x)] one has p'^irf) — pirf) € iM+. Therefore, it 
follows from the Cauchy-Riemann equations that for A in portions of Cl and Cr close enough (independently 
of Hm) to the interval [0, zyl(x)] one has 



(77) 



3?(t(A)) < 



for A on both Cl and Cr. Furthermore, it follows from the fact that p°(?/) G iR+ that 5R(— i0°(A)) < for 
A G Cl and 5R(i6lO(A)) < for A G Cr. Therefore 

(78) 3fi(T(A) - ^6l"(A)) < 
for A G Cl near the portion of Cm below iA{x), and 

(79) 3?(t(A) +i6i"(A)) < 

for A in the analogous portion of Cr. Next consider the situation to the left or right of the portion of Cm 
lying above the endpoint A — iA{x). Since 9{X) = and 5R(0(A)) < for A G [iA{x), iA\ we see that for Cl 
and Cr close enough (again independently of h^) to this part of Cm we again find that we have (|7|) on Cl 
and (^9|) on Cr. This shows that the jump matrix on both contours Cl and Cr is an exponentially small 
perturbation of the identity for small positive Hn, pointwise in A bounded away from the origin and lA. 
For A G [0, iA{x)], the jump matrix for N(A) factors (recall 0(A) = here): 

(80) 



exp 



i9{X)\ 



exp 



-I exp 



»g(A) \ 



-I exp 



i9{X) 



I'M 



Let Ll and Lr be two boundaries of a lens surrounding [0,i^]. See Figure ^. Using the factorization 
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Figure 4. Introduction of the lens boundaries and L^. 



we new define a new matrix function 0(A). In tlie region between and Cm set 



(81) 0(A) := N(A)cri 
In the region between Cm and Lji, set 

(82) 0(A) := n{X)a'~ 



1 i cxp 




1 



1 —i exp 







1 



Elsewhere in the upper half-plane set 0(A) := N(A). And in the lower half-plane define 0(A) by symmetry: 
0(A) =a20(A*)V2. 

These transformations imply jump conditions satisfied by 0(A) on the contours in Figure ^ since the jump 
conditions for N(A) are given. For A e L^, we have 

Hn J 



(83) 



0+(A) = 0_(A)ai 



1 —i exp 







1 



which is an exponentially small perturbation of the identity except near the endpoints. And for A € we 
have 



(84) 



0+(A) = 0_(A)ai^ 



1 —i cxp 
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which is also a jump that is exponentially close to the identity. For A e [0, iA{x)] we get 

" i 

(85) 0+(A)=0_(A) 

i 

as a consequence of the factorization (pO|). Since 0(A) := N(A) for all A in the upper half-plane outside the 
lens bounded by Ll and i^, we see that 0(A) satisfies the following jump condition on Cl- 



(86) 0+(A)=0_(A)ai 
the following jump relation on Cr: 

(87) 0+(A) = 0_(A)a^ 



1 







I exp 



I exp 



-(A) -^6l"(A) 



1 



and the following jump relation on the imaginary interval [i74(a;),iA] C Cm- 

1 



0+(A)=0_(A)ai 



I exp 



^(A) 



1 



All three of these matrices are exponentially close to the identity matrix pointwise in A for interior points of 
their respective contours. 

The matrix 0(A) is related to N(A) by explicit transformations. However, taking the pointwise limit of 
the jump matrix for 0(A) leads us to the following ad hoc model problem. 

Riemann-Hilbert Problem 3 (Outer model problem). Find a matrix 0(A) satisfying: 

1. Analyticity: 0(A) is analytic for A G C \ /, where I is the imaginary interval [—iA(x),iA{x)]. 

2. Boundary behavior: 0(A) assumes boundary values that are continuous except at X = ±iA{x), 
where at worst inverse fourth-root singularities are admitted. 

3. Jump condition: For X G I, 



(89) 



6+(A) = 6_(A) 



i 

1 



4. Normalization: 0(A) is normalized at infinity: 
(90) 6(A) ^ I as X^oo. 

This model problem is easily solved explicitly. 
Proposition 8. The unique solution of this Riemann-Hilbert problem is 

1 



(91) 



0(A) 



2i?(A)/3(A) 



R{X) - A - iA{x) R{X) + A + iA{x) 
R{X) + X + iA{x) R{X)-X-iAlx) 



where R{X)'^ = A^ + A{x)'^ and 
(92) 



PixY 



X + iA{x) 
X ~ iA{x) ' 



with both functions R{X) and /3(A) being analytic in C\I , normalized according to i?(A) 
as X oo. 



-X and /3(A) ~ 1 
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Using this matrix, we define an "outer" model for the matrix N(A) as follows. The idea is to recall the 
relationship between the matrix N(A) and 0(A), and simply substitute 0(A) for 0(A) in these formulae. 
For A in between Ll and Cm, we use (ISl]) to set 



(93) No,t(A) :-0(A)ai^ 
For A in between Cm and L;^, we use ( p2| ) to set 

(94) Nout(A) 6(A)a^^ 



-I exp 



tern. 

Hn J 



1 ■ ( 







1 



For all other A in the upper half-plane, set Nout(A) := 0(A), and in the lower half-plane set Nout(A) 
C2Nout(A*)*cr2. The important properties of this matrix are the following. 

Proposition 9. The matrix Nout(A) is analytic for all complex A except at the contours Ll, Lji, the imag- 
inary interval [0,iA{x)], and their complex-conjugates. It satisfies the following jump conditions: 



(95) 



Nout,+ (A) =Nout,-(A)(7i^ 



for Xe Ll; 



(96) 



Nout,+ (A) = Nout,-(A)ai 



(97) Nout,+ (A)=Nout,-(A)ai^ 



exp 



1 i exp 



^^(A) 



I'M 



ie{x) 



for Xe Lr, 



exp 



iejX) 



for Xe[0,iA{x)], 



with the jump matrices on the conjugate contours in the lower half-plane being obtained from these by the sj/m- 
metryN outiX*) = o'2Nout(A)*(T2. In particular, note that for X G [0,iA{x)], iwe /laue Nout,-(A)""'^Nout,+(A) = 
N_(A)^^N+(A). Also, if D is any given open set containing the endpoint A — iA{x), then Nout(A) is uni- 
formly bounded for A G C \ (13 U D*) with a bound that depends only on D and not on Hn . 



7. Local Analysis 

In justifying formally the local model Nout(A), we ignored the fact that the pointwise asymptotics for the 
jump matrices for 0(A) that we used to obtain the matrix 0(A) were not uniform near the origin or near 
the moving endpoint A = iA{x). We also neglected the breakdown of the asymptotics for aL{X), aii(A), 
and aM(A) near the points A = and A — iA. Consequently, we do not expect the outer model Nout(A) 
to be a good approximation to N(A) near A = 0, A = iA(x), or A = iA. In this section, we examine the 
neighborhoods of these three points in more detail, and we will obtain accurate local models for N(A) in the 
corresponding neighborhoods. 

7.1. Local analysis near A = 0. 



ASYMPTOTICS OF SEMICLASSICAL SOLITON ENSEMBLES: RIGOROUS JUSTIFICATION OF WKB 



21 



7.1.1. Local behavior of the matrix elements aL(A), qr^X), and aM(A). . Let e and S be small scales tied to 
such that fi-AT <C (5 <C e <C 1 as | 0. Let L be defined as the unique integer for which exactly N — L oi 
the numbers Aq, . . . , Xn-i lie strictly below ie on the positive imaginary axis. We want to compute uniform 
asymptotics for S{X) defined by (P) for XeCL^CR, and for T{X) defined by for A G Cm when | A| < <5. 



Lemma 4. When 5(A) > anc? |A| < 6 and with L defined as above, 

(98) exp (^-J2h(X)^=i + o(^^ 
Proof: We recall the integral formula (c/. equation (p8|)) 

(99) /fe(A) = -^ / dm dC d^giX,0, 
in which we expand the integrand in partial fractions: 



^'""^ ^(^'^)-A + e(0 ■ A-e(0 (A + e(e))^ ' (A - e(e))^ 

Since 5(A) > 0, for — hjy /2 < ^ < m^ + ftAf/2 and fc = 0, . . . , i — 1, we get 

(101) ^<^<.^< .. ' ...... 



|A + e(e)| - |A-e(e)| " K<^-e(0| " \i5 ~ e{mk ~ Kn \\m{5) - mk + Hn I2\ 
For such ^ we therefore have 

(102) 5(A,C) = 



|to((5) - mfe + ;iAr/2|2 



so summing over k gives 



(103) E4(A) = of;.^A.xi^x^ — ' , /9i2 l-o^^r ^^^1=^^^'" 

^ V ^ l"^((5)-mfc + ?l^r/2|■^y \^ (™ - "i((5))^ y 

because (S <C e, which proves the lemma. □ 

So only the fraction of terms /fc(A) with k > L contribute significantly to the sum for /(A). It is easy to 
check directly that exp(— /^(A)) is an analytic function for |A| < 5 whenever < fc < i — 1, so it makes no 
difference in these terms whether it is L'^{X) or X^,(A) that appears in the definition of Ik- Therefore, the 
terms in 5(A) and r(A) that can be significant for A near the origin are thus 

\k=L 



(104) si'Hx) := ( n ^-^ L (ie(™)(A) - i°-e(™)(A)) dm 



and 

,105, rriA):. (n (^f *''(i:,„,(A,-I».„,A)) .2co» (^) ^ 

Here we have written the integrals in the exponent using the change of variables m = m{rj). So Lemma ^ 
simply says that S'(A) = sf\x){l + 0{hi^ I e)) and T(A) = T[°\x){l + 0{hN / e)) uniformly for |A| < 5. When 
A is close to the origin along with the points Afc contributing to T(A), the ladder of discrete nodes appears 
to become equally spaced. The next lemma shows that this is indeed the case. 

Lemma 5. Let X^-k for A; — 1, 2, 3, . . . be the sequence of numbers defined by the relation: 
(106) AA._,.:^--^(fc~l/2), 
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which results from expanding the Bohr-Sommerfeld relation ( |Jq ^ for \N-k small, and keeping only the dom- 
inant terms. Define 

(107) 5f (A) := in ^-^j exp (i°'(o)™(A) - i"-e.(o)™(A)) dmj , 



and 



(108) 



,/c=L 

X 2 COS f^(^A-A) 
T/ien, /or 3(A) > and |A| < 5, 

(109) rf)(A) = (A) (i + o (^1^ log (^^) 

where we suppose that the scale e is further constrained so that the relative error is asymptotically small. If 
A is additionally bounded outside of some sector containing the positive imaginary axis, then 

,2 



(110) sr\X) = s'^i'\\)\\ + o\j- 

Proof: We begin by observing that ioi k — L, . . . , N — 1, the distance between \k and A^ is much smaller 

1 /2 

than the distance between Afe and Afe+i, as long as e <C hj^ . More precisely, we have 

(111) \\k - Xk\ ^ 0{h%{N ~ kf) . 
Decompose the quotients as follows: 

(112) = DiX)CiX)LiX) and = DiX)LiX) , 
where 

k=L ^ 



(114) 



C{X) :=cos (^^j'\\j^)d^ sec (^-^N - ^p\Q)X^ , 



/ 1 n7nL-\-hN /2 

(115) L(A) exp ( — ([leM (A) - l!!,(o)™(A)] - [1° eM(A) - 1° e'(o)™(A)]) 



dm 



and 



(116) L{X) exp ( — ([L°(™)(A) - L°,(o)™(A)] - [i",(„)(A) - L° ,,(o)„(A)]) rfmj . 

First we deal with L{X) and L{X). Since e(m) is difFerentiable and m is small we have e(m)— e'(0)m = 0(e)- 
Also, the interval of integration is 0(e) in length. Although the integrands in (115) and (11€) arc not 
pointwise small, upon integration it follows that 



eM , , „/e2 



(117) L{X) = 1 + O — and L(A) = 1 + O . . 

\hN J \hN 

uniformly for all A in the upper half-plane satisfying |A| < (5. Here we are assuming that e <C ft^ . 

For the moment, let's drop the conditions 3(A) > and |A| < (5 and instead consider A to lie on the sides 
of the square centered at the origin, one of whose sides is parallel to the real axis and intersects the positive 
imaginary axis halfway between the points A = Al and A ~ Al-i. Note that the estimate (111) implies that 
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the sides of the square intersect the real and imaginary axes a distance from the origin that is approximately 
e. Therefore the square asymptotically contains the closed disk |A| < 5 because 5 <^ e. We will show that 
for A on the four sides of the square, both D{X) and C(A) are very close to one. We write D{X) in the form 

Afc — Afc \ 



N~l 

(118) D{X) = 

k=L 




A- Afc 



First consider the top of the square: for 3(A) = —i{\L + Al-i)/2, we easily see that 
(119) |A-A.|>||,*-L + l/2) and ^ « (7 



for k = L, . . . , N — 1. Combining this with (111), we get 

(120) K^^o(^%i^->^n h^^of ^-(^-'^^^ 



A -A* V e / A-Afe \hNik- L + 1/2), 

Summing these estimates over k (it is convenient to approximate sums by integrals in doing so) , we find that 

Consequently, for A on the top of the square. 



IN \ "TV 



e 



(122) D{X)^l + 0[—log[ — 

\hN \hN 

An estimate of the same form holds when A is on the bottom of the square, where 3(A) = i(AL + Ai_i)/2. 
When A is on the left or right side of the square, so that |3?(A)| = — j(Ai, + Al_i)/2, both |A — A^,|^^ and 
|A — Afcl^-^ are 0{e~^). By the same arguments as above, we then have for such A that 

(123) D{X) = l + o(-^ 

\riN 

Now we look at C(A) on the same square. Generally, for such A which are of order e in magnitude, we have 

(124) C(A) = 1 + sec ( -nN - ^p°{0)x] . 
When A is on the top or bottom of the square, we have 



(125) 



sec ^-vrA^- ■^p''(0)A 



<1, 



and when A is on the left or right sides of the square, the same quantity is exponentially small. It follows 
easily that for A on any of the sides of the square, 

(126) C{X) = l + o(-^ 



So uniformly on the four sides of the square, we have 



(127) D(A)C'(A) = 1 + O — log 



---o ( ^ ) 

But the product Z3(A)C(A) is analytic within the square, so by the maximum principle it follows that the 
same estimate holds for all A on the interior of the square, and in particular for all A in the upper half-plane 
with |A| < S. This shows that 

(128) rf)(A) = Tf (A) (1 + O (^1^ log (^^ 

holds for all such A. 
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Now to control the relationship between S'^"-' (A) and Slf' (A) we consider A to lie outside of some symmet- 
rical sector about the positive imaginary axis, of arbitrarily small nonzero opening angle 2a independent of 
Hn. Since 3(A) > 0, we get 

(129) |A-A;.|>|A-A..|>J^ = '^"y„-^7'/^' . 

sm(Q;) p"(0)|sm(a)| 



Combining this result with (111), we find 

(130) ^0{hN{N-k)) and ^ OihN{N ~ k)) . 
A — A^. A — Afc 

Summing these estimates over k one finds that 

(131) D{\)^l + o(f- 



Combining this with the estimate (117) of L(A) — 1, we find that 



(132) ^(o)(A)^4")(A)(^l + 0(^— 

for all A in the upper half-plane with |A| < 5 and bounded outside of the sector of opening angle 2a about 
the positive imaginary axis. This completes the proof. □ 

Without any approximation, s!f^ (A) can be rewritten in the form: 



(133) (A) 



(-tC)-'^ aC) r(l/2 + zQiN + zQ^+^'^TiN + 1/2 - zQ 
r{l/2~iO{N-iC)^~<T{N + l/2 + i(:) 

andTf)(A) can be rewritten in the form: 

(134) r(°VA) - ( ,. aiv + <)^+<r(7v + i/2-<) 

where N := N — L and we are introducing a transformation (po to a local variable ( given by 

(135) C = ipoW-=-^P°iO)X/hN- 

These formulae come from evaluating the logarithmic integrals exactly, which is possible because e(m) has 
been replaced by the linear function e'(0)r7i, taking advantage of the equal spacing of the A^^ to write the 
product explicitly in terms of gamma functions, and then using the reflection identity for the gamma function 
to eliminate the cosine from T2^\x). Now, the integer N is large, approximately of size e/fiN. But for |A| < (5, 
N is asymptotically large compared to C because S e. These observations allow us to apply Stir ling- type 
asymptotics to iSj^"* (A) and Tj^-* (A) . 

Lemma 6. In addition to all prior hypotheses, suppose that S'^ <C sHn- Then, 
and 



(136) ^r(A)-e2<(-zC)-<(zC)-^«-)-^-— l + O, . 

r(i/2 + ig V X^^N 



r(i/2-iC)2 V V^N, 

Proof: Asymptotically expanding the gamma functions for large N , we find that 

(138) 5f (A) = e^<(-zC)-^«)--^ ra/2-<) ' ^^^'"^^ ' (' + (^)) 

and 
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where 

140 A(C,A^) := — — = — -^^ 

^ ' ^ ' (7V + <+ l/2)^+< [N-iCY-< 

Next, expanding A(C, iV), one gets worse error terms: 



(141) A(C,7V) 




5 V 

N 



Combining these estimates and noting that 1/N — 0{hN /e) completes the proof of the lemma. □ 
With these results in hand, we can easily establish the following. 

Proposition 10. Let A be in the upper half-plane, with |A| < ft^, where 3/4 < a < 1, and let A be bounded 
outside of some fixed symmetrical sector containing the positive imaginary axis. Then 

(142) SiX) = e^'C(_,^)-<(,^)-< r(l/2 + g ^ o(4"/3-i)) , 

where ^ — tpoW ■— ~ip^{0)X/hN. 

Proof: According to Lemmas |[ ^, and ^, the total relative error is a sum of three terms: 



(143) O I — ) and O ( ^ ) and o' ^ 



2 , 



Note that since Hn <C (5, the order /f- term is always dominated asymptotically by the order 5 /eh^ term. 
The error is optimized by picking e so that the two possibly dominant terms are in balance. This forces us 
to choose e ~ S'^^'^ . The proposition follows upon taking 5 — h%. □ 

Proposition 11. Let A be in the upper half-plane, with |A| < h^, where 3/4 < a < 1. Then for all v > Q, 
however small, 

(144) r(A) = V + ^(^^ n ' 
where C, — (poiX) := —ip'^{0)X/hN . 

Proof: In this case, according to Lemmas ||, ^, and H, the total relative error is a sum of three different 
terms: 

(145) o(^) and o(^\og(^)) and O ^ 



Again, since Hn <C S, the order h^/e term is always dominated asymptotically by the order S'^/eh^ term. 
For any cr > 0, we have 

riN \riN J \nN Km 
so we can eliminate the logarithm at the expense of a sli ghtly larger error. Taking 5 = hf^ as in the statement 
of the proposition, and using the cruder estimate (146), the nearly optimal value of e to minimize the total 
relative error is achieved by a dominant balance between the right-hand side of (146) and the term of order 
5^/efi.Ar. The balance gives e = /i^, with 

(147) /5=^- 

3 + (T 

With this choice of e, the total relative error is of the order ft^, with 

, , 4a + 2(a - l)cr - 3 4 

148 2a- 1-/3= <^a-l 

3 + (T 3 

with the inequality following because cr > and a < 1. The inequality fails in the limit cr ^ 0. Therefore, 
for each arbitrarily small > 0, we can find a cr > sufficiently small that ^>Aa/'i — 1 — v. This gives us 
a slightly less optimal estimate of the relative error: simply 0{h'^^^ ^ '^), which completes the proof. □ 
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7.1.2. The model Riernann-Hilhert problem. To repair the flaw in our model Noui(A) for the matrix N(A) 
related to the nonuniformity of the approximation of the jump matrices near the origin, we need to provide 
a different approximation of N(A) that will be valid when |A| < h'^ for some a € (3/4, 1). The local failure 
of the approximation is gauged by the deviation of the matrix quotient N(A)Noiii(-^)^^ from the identity 
matrix near the origin. It turns out to be more convenient to study a conjugated form of this matrix (which 
also deviates from the identity for A near the origin). Namely, for |A| <K^, set 

(149) F(A) := e-'^(°)'^=*/(2ft^)ai^ (ic7i)6(A)-iN(A)Nout(A)-^6(A)(-iCTi)£7i^e'^W'^^/(2ft^' 
if K(A) < and 

(150) F(A) := e-'^W'^^/(2'^~Vi^6(A)-iN(A)Nout(A)-^6(A)ai^e*^W'^='/(2ft;.) 

if 3?(A) > 0. It is easy to check that as a consequence of the boundary conditions satisfied by the matrix 0(A) 
on the imaginary axis near the origin, the conjugating factors arc analytic throughout the disk |A| < K^. 

For later convenience, let us assume without loss of generality that the auxiliary contours Cl, Cr, Ll 
and Lij are straight rays in some /lAr-independent neighborhood of the origin. It is easy to write down the 
jump conditions satisfied by F(A) on these four rays and also on the positive imaginary axis. We find: 



(151) 



F+(A)=F_(A) 



1 aL(A)e-^^W/'^« 
1 



for A e Ci 



(152) 



F+(A)=F_(A) 



1 

aii(A)e^^(°)/'*" 1 



, for A e Cfl , 



(153) 



F+(A)=F_(A) 



1 

.ie-'(.f){x)-e(o))/h^ ^ 



, for Xg Ll . 



(154) 

and for A e Cm, 
(155) 

F+(A)=F_(A) 



F+(A)=F_(A) 



1 -ieiioW-om/hN 
1 



for A e L 



R ■ 



l + [mM(A)+e^(^)/'^"] 



1- [iaM(A) + e'^(^)/'^"] 



The jump relations satisfied by F(A) on the complex conjugate contours in the lower half-plane follow from 
these by the symmetry F(A) = cr2F(A*)*o-2. 
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Now, for A e Cl with |A| < h 



AT; 



4>{X) + i{9{X) - 6(0)) - 2i0"{X) 



I exp 



■N 



(156) ^(^(A) - e{0)) - 2i(0"(A) - 6*0(0)) 

= zexp 

= te^^SiX) {l + O {h%^-')) 




^-^'^(-<)-'^«)-^ ra/2-S + 



"Xi/2-*C) 

where m the last hne ( — fo{X) with the change of coordinate being given by (135). In these steps, we 
used the relation (75), the fact that from Proposition || we get ^(A) = 0, and, according to ( |ll| ) and the 
quantization condition ( |T^ ) on Hn, 2e°{0)/hN = 2TrN e 2tt1. We have also used the fact that 

/ , 1 f d9° , , d9 , ,\ 



which follows directly from the definition ( ^A\ ) of 6*" (A), the definition ( |69| ) of 0{X), and the relation (|65|). In 
a similar way, for A G Cr with |A| < we get 

(158) a«(A)e^'^(")/'-- = ze(^''-)'^(-zC)--'^(»C)-''^ r||^' ^ + ^ (4"^'"')) ' 
and for A G U Ca/ U with |A| < h%, 

(159) g±.(e(A)-fl(o))A„ ^ gT-C (1 + o , 

with C = ipo{X). Finally, when A G Cm and |A| < ft^ we have for arbitrarily small v > 0, 
mAf(A) + e'^(^)/'"'" = [1 - r(A)] 



= l-T(A) 

(160) _ _ 27re2<(-<r;2iC 

^'^{l/2~J0^ 



^(l/2-^C)2 ' ^ V ^ 

again with C = (y5o('^)- The last step follows because 27re^*^(— iC)^^*''/r(l/2 — iQ"^ is uniformly bounded on 
Cm- 

Let Cl, Ci{, Li, Lfl, and Ca/ denote the straight rays that agree with the corresponding contours in a 
fixed neighborhood of the origin in the A-plane, but lying in the C-plane (according to (135), C is a simple 
rescaling of A by a positive number). These rays are oriented contours, with the same orientation as the 
original contours: Cm, Ll, and Lji are oriented outwards from the origin toward infinity, and Cl and Cr 
are oriented inwards from infinity toward the origin. Let the union of these contours with their complex 
conjugates be denoted So- Consider the following Riemann-Hilbert problem. 

Riemann-Hilbert Problem 4 (Local Model for the Origin). Find a matrix F{C) with the following prop- 
erties: 

1. Analyticity: F{C) is analytic for C G C \ Sq. 

2. Boundary behavior: F{() assumes continuous boundary values on Eq. 
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3. Jump conditions: The boundary values taken on Eo satisfy 
(161) 



(162) 



(163) 



(164) 



F+(C)=F_(C) 



F+(C)=F_(C) 



^ ^ ^ r(i/2-zC) 







^ ^> r(i/2-iC) 



F+(C) = F_(C) 



1 



F+(C)=F_(C) 



and 



(165) F+(C)=F_(C) 



1 -ie-""^ 
1 

27re^<(-zC)-^'^ 

r(i/2~*C)2 



2TTe^<{-iO-'^"^ _ ■ 

r(i/2-*c)2 . 



le 



r(i/2-*c)2 

27re'<(-<)-^< 

r(i/2-iC)2 



C e Cm . 



On t/ie contours in the lower half-plane, the jump conditions are implied by the symmetry F(^*) = 
a2F(C)V2. 

4. Normalization: F(C) is normalized at infinity: 
(166) F(C) I as C ^ oo • 

Unfortunately, we cannot solve Riemann-Hilbert Problem ^ explicitly. Luckily, we will not require an 
explicit solution. However, existence of a solution is an issue that must be resolved, and we need to obtain 
a decay estimate that quantifies the normalization condition (^66|). These questions are addressed via the 
abstract theory of Riemann-Hilbert problems. 

Proposition 12. The local model Riemann-Hilbert Problem ^ has a unique solution satisfying F{() = I + 
0{1/C) as ^ — > oo, uniformly with respect to direction. Also, det(F(C)) = 1. 

Proof: Each Riemann-Hilbert problem is equivalent to an inhomogeneous system of linear singular integral 
equations. It must be shown that the matrix singular integral operator associated with these equations is of 
Fredholm type, with index zero. Then it must be shown that there are no homogeneous solutions, at which 
point one has existence and uniqueness of a solution to the inhomogeneous system. Finally, one maps the 
solution of the integral equations to the unique solution of the Riemann-Hilbert problem and it remains to 
verify the rate of decay to the identity matrix as C oo. 

The theory we will use is the theory of matrix Riemann-Hilbert problems on self-intersecting contours, 
with boundary values taken in spaces of Holder continuous functions. This theory is summarized in a 



self-contained way in the appendix of |KMMOO|. 

The first step is to establish that the operator of the associated system of singular integral equations is 
Fredholm index zero on an appropriate space of functions. As described in | KMMO(| , this follows from two 
facts. First, on each ray of the contour Eg the jump matrix Vp(C) := F_(C)~^F+(C) is uniformly Lipschitz 
with respect to (, and differs from the identity by a quantity that is 0{1/Q for large Second, the limiting 
values of the jump matrix, taken as C — > along each ray of Eg, are consistent with a bounded solution F{() 
near C = 0. This means the following. Suppose that F(C) has a limiting value, say a matrix Fq, as C ~* 
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inside one of the sectors of C \ Sq. Using the hmiting value of the jump matrix Vp(C) at the origin along 
one of the rays of Eq bounding that sector, one can compute the limiting value of F(C) at the origin in the 
neighboring sector. This procedure can be continued, moving from sector to sector of C \ Sq in the same 
direction, until one arrives once again in the original sector, with a matrix Fi. The consistency condition is 
simply that Fi = Fq, which upon elimination of Fq can be viewed as a cyclic relation among the limiting 
values of the jump matrix Vp(C) taken along each ray of Eg as C — > 0. It is easily checked that this cyclic 
relation indeed holds for Riemann-Hilbert Problem ^. 

The second step is to establish existence and uniqueness of the solution ¥{()■ The fact that the associated 
singular integral equations are Fredholm index zero means that, in a certain precise sense, the Fredholm 
alternative applies to our Riemann-Hilbert problem. The inhomogeneity is the normalization to the identity 
matrix at C = oo. The corresponding homogeneous Riemann-Hilbert problem has exactly the same form 
except that the normalization condition is replaced by the condition F(C) — > as C ^ oo. We will have a 
unique solution of Riemann-Hilbert Problem ^ if it can be shown that no such homogeneous solutions exist. 
For this purpose, it is sufficient that the jump matrix Vp(C) should have a certain symmetry with respect 
to Schwartz reflection through the real axis in the ^-plane. For the orientation of Eq described above, the 
required relation is: 

(167) Vp(C)-i=Vp(C)t 

for ( e Eo n C+. It is easily checked, using the symmetry Vp(^*) — <T2'Vp{C)*a2 and structural details of 
Vp(C) in the upper half-plane, that this relation holds, and this means that Riemann-Hilbert Problem ^ has 
a unique solution F(C). 

The third step is to establish that the unique solution F{C) decays to the identity for large C like 1/C- 
The Holder theory that we have been using generally provides a solution F(C) under these circumstances 
that takes boundary values on Eq that are Holder continuous with exponent fi and that differs from the 
identity matrix by 0(1/C^) as C — > oo, for all /i strictly less than 1. This fact can be traced to the compact 
embedding of each Holder space into all Holder spaces with strictly smaller exponents. The compactness 
is needed to establish the Fredholm property of the Riemann-Hilbert problem. So to obtain the required 
decay, we need an additional argument. The condition that is required to obtain the 0{l/() decay is that a 
signed sum of the mean values of ( ■ (vp(() — I) taken as C ^ oo along each ray of Eq (the signs are related 
to the orientation of the individual rays) is zero ||KMMOdt . Now along each ray of Eq except for Cm and its 
conjugate {i.e. the imaginary axis in the ('-plane), Vp(C) decays to the identity exponentially fast as C ^ oo. 
So these rays to do not contribute to the sum and it is only necessary to check the imaginary axis. When 
C e Cm, 



(168) C-(vp(C)-]I) = ^ 



as C — > oo, and for ( on the negative imaginary axis oriented upwards, 

-1 -ie""^ 



(169) C-(vf(C)-]I) = 3^ 

as C — > oo. The limits of these quantities do not exist as C ^ oo due to the oscillations on the off-diagonal. 
But the mean values exist and are equal, and it turns out that they enter the sum with opposite signs due 
to the orientation of the contour rays. Thus, the required sum of signed mean values indeed vanishes. This, 
along with the analyticity of the jump matrix Vp(C) along each ray of Eg establishes that F{C) — I = 0(1/C) 
as C — > oo. 

Finally, we check that det(F(C)) = 1. Taking determinants in the jump relations we see that on all rays 
of the contour, det(F+(()) = det(F_(()). Since the boundary values taken by F{C) on Eq are continuous, 
we discover that det(F(C)) is an entire function. Since this function tends to one at infinity, it follows from 
Liouville's theorem that det(F(C)) = 1. This completes the proof of the proposition. □ 



7.1.3. The local model for N(A) near A ~ 0. From (149) and (|150| ) we can express N(A) in terms of F(A) 
for |A| < /i^. For 5ft(A) < we have 

(170) N(A) = 6(A)(z(Ti)ai^e*^(o)"3/(2'i")F(A)e-^^(°)"^/(2;i«)^ V (_iai)6(A)-iNout(A) 
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and for 5i(A) > we have 

(171) N(A) = 6(A)ai^e'^(°)'^^/(2'^^)F(A)e-^^(°)'^=*/(2ft«)^T^Q(;^^-iN^^^;^) _ 

To obtain a local model for N(A) near the origin, wc simply replace F(A) in these formulae by the approxi- 
mation F{tpo{X)). With |A| < h%, we set for 3fJ(A) < 0, 

(172) Norigin(A) := 6(A)(i<7i)ai^e'^(o)-^/(2ft«)F(^o(A))e-'^(o)-^/(2ft«)(7i^(-ic7i)6(A)-iNo„t(A) 
and for K(A) > we set 

(173) NoHgi„(A) := 6(A)<7i'^e'^W<^='/(2'^^)F(^o(A))e-'^(o)'^='/(2'^^)ai'^6(A)-iNout(A) . 
The most important properties of this matrix function are easily seen to be the following. 

Proposition 13. The m,atrix Norigin(A) is a piecewise analytic function of A in the disk |A| < ft^, with 
jumps only on the locally straight-line contours Cl, Cr, Ll, Lr, and Cm, and their conjugates in the lower 
half-disk. The jump relations satisfied by Norigin(A) on these contours are the following: 

(174) 

1 



^origin,— (^) NQj-igin,+ (A) = (7i 

(175) 

^origin,— ('^) Norigin,+ 

(A) = cr^ ^ 







e I *g i«g r(i/2-iC) 



(176) 

(177) 

(178) 
where 

(179) 

with 
(180) 



-(A) ^Norigin,+ (A) = CTi 



-(A)-'Norigin,+ (A) =ai 



1 i[e-iS{X)/hN _ g7rfg-ie(0)/ftN-) 



I ^^gi6l(A)/Sjv _ g-7rCgi6l(0)/;tjv-) 



1- J 
^1= , 



'1 5 



XeLn, 



^0 1 

Norigm,-(A)"^Norigin,+ (A) = " v(A)cri " , A € Cm , 

gie(A)/ftjv(i + (1 - e-'^«-'(^(^)-^(o))/^")Z) , 
iZ(e''^+'(^(^)-^W)/'*^ + e-^C-'(«W-^(o))/ftJv _ 2) , 
i + iZ, 

g-ie(A)/ft„^l + (1 - e''^+'(^(^'-^W)/'^'^)Z) , 



^^21 
^^22 



r(l/2-zC)2 

anrf where ( = (po{X). The jumps on the corresponding contours in the lower half-plane are obtained from 
the symmetry Norigin(A*) = cr2Norigin(A)*cr2. The matrix Norigin(A) is uniformly bounded for \X\ < K^, with 
a bound that is independent of h^- Also, det(Norigin(A)) = 1 and when \X\ = h'^, 

(181) Norigin(A)No„t(A)-l = I + 0(/i]v"") • 
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7.2. Local analysis near A = iA. 

7.2.1. Local behavior of ai(A), aj?(A), and qmW- As before, we suppose that e and S are small scales 
satisfying fij^ <C ^ <C e ^ 1 as Hn tends to zero. We redefine the integer L so that exactly the first L of the 
numbers Aq, . . . , Xn-i lie on the positive imaginary axis above i{A — e). We will suppose that 3(A) < A, 
and |A — iA\ < S and we will deduce asymptotic formulae for T(A) given by (|3^) valid for such A, and for 
5(A) given by ( pl| ) when A is also bounded outside of some downward-opening sector with vertex at iA, in 
the semiclassical limit Hn 0. First, we establish a result that is the analogue of Lemma ^. 

Lemma 7. When 3(A) < A and |A — iA\ < S and with L defined as above, 
(182) exp 



(-g4.))...o(^ 



Proof: We again estimate /^(A) using the integral formula ( p9[ ) with integrand g{X,^) given by (100). 
Given our conditions on A, for mk — hN/2 < ^ < mk + /iAr/2 and k > L we have 



(183) 



11 1 1 ^ / 1 

< 7^ 7TT7 < T-^. 7TTT < VT^. ^ Z TITTT = O 



|A + e(0| - |A - e(0| " HA - S) ^ e(e)| " \i{A -5)^ e{mk + hN/2)\ \ \m{iA - iS) - - 
For all such ^ we therefore have the estimate 

(184) 5(A,0=O^ ^ 



|m(zA — iS) — mk — ^jv/2| 



Summing over k gives 

/ ^-1 2. \ f rm(iA-t) \ 

(185) E ^'^(A) - o [nl g 1^^^^ _ _ _ j = O (^ft. M^X^TjpH^j 

which is 0{his! /e) because 5 ^ e, and the lemma is proved. □ 

As was the case when A was near the origin, only certain terms are important when A — zv4 is small, as 
a direct consequence of Lemma ^ and the fact that exp(— /jt(A)) is analytic for such A when k > L. The 
important terms when |A — iA\ < 6 are 

(186) S['^\X) := ( n ^-^] exp C U'^^^^W ' i° e(™)(A)) dm] , 



and 



(-) -rw - (n (;^/l„,(^."..<^>-^-.<^))-] --K^) 



Lemma states that S{X) = S'f^^(A)(l + 0{hN/e)) and r(A) = Tf^^(A)(l + 0(ftAr/e)) as ^iat ^ for 
\X-iA\ < S. 

The analogue of Lemma |^ says that for X — iA small the sequence of numbers Xq, . . . , Xj^^i contributing 
to ^(A) and T(A) can be replaced essentially by a "straightened-out" sequence with uniform density. 

Lemma 8. Let Xk for fc = 0, 1, 2, . . . be the sequence of numbers defined by the relation: 
(188) Afe:=zA+-j^(fc + l/2), 
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which results from expanding the Bohr-Sommerfeld relation for Afe near iA and keeping only the domi- 
nant terms. Define 



(189) 



and 



k=o 

X exp 



fe=0 ^"^^ 



T^'w ■■= n 



(190) 



I'M 



X exp I -5- 



(^LiA+e'{M){m-M)W ~ ^-jA~e' (M)(m-Af) ('^)) "^"^ 



T/ien, for 3(A) < A anrf |A - iA| < J, 

(191) ^)(A) = Ti'^)(A) (l + O log 

where we suppose that the scale e is further constrained so that the relative error is asymptotically small. If 
A is additionally bounded outside of some downward opening sector with vertex at iA, then 

(192) S['^\X) = 5^^^^ (A) (^1 + O (^ii 

Proof: The proof of this Lemma follows that of Lemma ^ almost exactly and will not be repeated here. 
The only difference is that the square in that proof should be replaced here by the rectangle whose top side 
is 3(A) = A and — e < 5R(A) < e and whose bottom is 5(A) = —i{XL-i + Al)/2. □ 

Without any approximation, S2^\X) and T2"^\x) can be rewritten in a more transparent form by ex- 
pressing the products in terms of gamma functions and evaluating the logarithmic integrals exactly. Introduce 
a local variable ^ in terms of a transformation ipiA given by the relation 

(193) ^ = ^^^(A) :=p"(zA)^-^, 
and let B be the positive constant 

, , 2iAp"(iA) 

(194) B := — 

In terms of these quantities, one finds that 5*2*"^ ■'(A) and T2*'^''(A) take a simple form: 

(195) S^^^\X) = (-*C)^fr(l/2 - *C) • F(C, B, L) ■ W{C. B, L) , 
and 

.(iA)i^^ _ MKY ^ 
r(i/2 + zC) 

where 

r(B-< + i/2) 



(196) n^^^ (A) - ■V{C,B,L)-W(C,B,L)., 



(197) V{CB,L) 



T{B - L-i(+ l/2)r(L -iC+ 1/2) 
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and 

Now, B '3> L '3> \C\ because B is proportional to and L is of the order of e/h while |(^| — O^S/Hn). 
So again we can use Stirling's formula to extract the dominant asymptotic contributions to iS'2*'^''(A) and 
T2^^\X) as Hn tends to zero. 

Lemma 9. As Hn tends to zero through positive values, 



(199) S^'^\X) = ^e-<(-*C)-^^r(l/2 -^Q).(l + 



'2 

and 



ehN 



.HA),,, V2^e-<(zC)< A , 5 



^2 



Proof: Using Stirling's formula, one expands V{C,, B, L) to find 

(201) ViC S, L) - ^ _ i _ + l/2)^-^-C(L - zC + l/2)^-C + ^ (, V 

The error here is dominated by the fact L ~ e/^Ar is the smallest large number involved. Now we expand 
the powers that remain in conjunction with those in VF(C, B, L) to find 

(202) nC, DWiC i., L) . ^ (l + O (^) ) (l + O (|- 



Since S <^ Hn, the relative error is dominated by 0{6^ /ehiy). Using this expression for the product VW in 
(195) and ( |196| ), the lemma is proved. □. 

In exactly the same way as in our study of the local behavior near the origin, we may combine Lemmas 0, 
^, and ^ and choose the "internal" scale e in terms of S and Hn to obtain asymptotics for S{X) and T{X) 
with optimized relative error. 

Proposition 14. Let 3(A) < A, with \X — iA\ < h'^, where 3/4 < a < 1, and let X he hounded outside some 
fixed symmetrical sector with vertex at iA and opening downward. Then 

(203) 5(A) = ^e-<(-zC)-^^r(l/2 - C) (l + O (h'^'^-')) , 

where C, — ifiA{X). 



Proposition 15. Let 3(A) < A, with \X — iA\ < h'^, where 3/4 < a < 1. Then for all > 0, however 
small, 

where C = ipiA{X). 

7.2.2. Why a local model near X — iA is not necessary. In particular it follows from these considerations that 
both S{X) and T(A) are uniformly bounded functions on their respective contours in any fixed neighborhood 
U oi X ~ iA. We claim that the quotient of the jump matrices for N(A) and Nout(A) is uniformly close to 
the identity matrix in U as Hn — > 0. Since Nout(A) is, by definition, analytic throughout U, it suffices to 
show that the jump matrix vn(A) := N_(A)^^N-|_(A) is uniformly close to the identity in U. This will be 
the case if ai(A), afl(A), and ajv/(A) are uniformly small on their respective contours. Now for A G C'm, 

(205) aMiX)=iexp I^]t{X), 
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SO, with 0(A) being real and strictly negative for A e Cm nU according to Proposition |^ and T(A) being 
bounded, we see that cm (A) is in fact exponentially small as Hn tends to zero through positive values. 
Similarly, for A G C^, 

(206) a.(A) = . exp (^^M^) S(A) = . =p f^^^^^^) S,X) . 

where we have used ( [ts] ) and the fact that 9{X) = on Cm above A ~ iA{x). Since d'^{iA) = 0, it is 
possible to choose the neighborhood U small enough (independent of Hn) so that 5R((^(A) — 2i0°{X)) < 
throughout U. Since S{X) is bounded, it then follows that for A € Cl n C/, aL(A) is exponentially small as 
Hn tends to zero through positive values. Virtually the same argument using (^ in place of shows 
that afl(A) is also exponentially small for A e Cr D U (it may be necessary to make U slightly smaller to 
have 5R(0(A) + 2i6^{X)) < throughout U). It follows that vn(A) — I is exponentially small uniformly for 
the contours within U. 

For this reason we expect that the outer model Nout(A) will be a good approximation to N(A) near X — iA 
even though S{X) — 1 and T(A) — 1 are not small. We do not need to construct a special-purpose local model 
for N(A) in this case. 



7.3. Local analysis near A = iA(x). Let Z3 be a circular disk centered at X = iA of sufhcicntly small 
radius (independent of Hn) that ^ D and iA ^ D, and that Ll and each have exactly one intersection 
with dD (of course Cm will have two intersection points with dD). This situation is possible as long as 
X ^ 0. The case of a; = is a degenerate case that we will not treat in detail here. 

Since x ^ and therefore A{x) < A, it follows from the definition (^) of p{ri) that as A tends to iA{x) 
along /, 6{X) vanishes like (A — iA{x))^^'^, and not to higher order. Since 6{X) may be extended from / to 
be an analytic function in D except for a branch cut along the part of Cm in D lying above the center, and 
since p(A) extended to this cut domain from / is nonzero, it is easy to see that the function 



(207) 



defines an invertible conformal mapping of all of D to its image. Consider the local variable C defined by the 
relation C = ^PiA(x){^)- The image of D in the ^-plane is a neighborhood of = that scales with hjs! such 

2/3 

that it contains the disk centered at C = with radius Ch^ for some constant C > 0. The transformation 



(207) maps /n D to a ray segment of the positive real ^-axis, and takes the portion of Cm in D lying above 
A = iA{x) to a ray segment of the negative real ^-axis. We suppose that the contours Ll and Lr have been 
chosen so that (piA(x)iLL H D) and <PiA(x){LR n D) are straight ray segments with angles — 7r/3 and tt/B 
respectively. 

For C G 'PiA{x){D), the matrix S(C) := ^ ^{'^iA(x)^0)'^i ^ satisfies the following jump relations: 



(208) 



s+(C) = s_(c) 



i 

1 



for C e M+ n LpiA{x){D), oriented from right to left, 
(209) S+(C) = S_(C) 



,3/2 



1 -ie-< 
1 

on the part of the ray arg(C) = — 7r/3 in ^PiA(x){D), oriented toward the origin, 

1 — le 



(210) 



s+(c) = s_(c) 







1 
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on the part of the ray arg(C) — tt/3 m (^j^j-^.) (I?), oriented toward the origin, and finally 

1 

(211) S+(C) = S_(C)' 



for C e M- n (piA{x)(D), oriented from right to left. The jump relation on the negative real ^-axis follows 
from the formula ( |6^ ) which applies because 0(A) — N(A) here. In ( |62|) one uses the fact that ^(A) = for 
A e Cm above iA{x), and the relation (|7^ ) giving 0(A) above iA{x) in terms of the analytic continuation of 
0{X) from /, which one writes in terms of the local coordinate ^. 

As ^ oo on all of the rays except for the jump matrix for S(C) decays exponentially to the identity 

matrix. These jump conditions were precisely the ones that were neglected in obtaining the outer model. 

i-j _ 1- J 

That is, the matrix S(C) :— cr^^ 0{(p~2f^,^^{C))o'i^ defined for C S ^iA(x){D) is analytic except on the 
positive real ^-axis, where it satisfies 



(212) 



s+(C) = s_(C) 



i 

1 



It follows that the matrix S(C) can be decomposed into a product of a holomorphic prefactor depending on 
h]^ and a universal {i.e. independent of hpf) local factor that takes care of the jump. We therefore may write 



(213) 



s(c) s'^°'(C)s'°'=(c) 



where S (C) is holomorphic in ipiA{x)iD) and where 



(214) 



s>-(c) -^i-cr^' 



' 1 1 ' 


1 


-1 1 


"71 



(_C)V4 (_C)l/4 



Note that S''°'(C) has determinant one. Its matrix elements are of size 0{hj^^^^) for ( e ipiji{x){D). It is 
easy to write down an explicit formula for S^°'(C) because both S(^) and S'°'^(C) are known. 
We will now approximate S(C) by 



(215) 



S(C) := S'>°'(C)S'°^(C) 



be the contour shown in 



where S (C) is the solution of the following Riemann-Hilbert problem. Let ^^^(j:) 
Figure ||. 

Riemann-Hilbert Problem 5 (Local Model for iA{x)). Find a matrix S^°'^{() with the following proper- 
ties: 

1. Analyticity: S'°'^(C) is analytic for ( \ Y,iA(x)- 

2. Boundary behavior: S'°'^(C) assumes continuous boundary values on T,iA(x)- 

3. Jump conditions: The boundary values taken on Tha{x) satisfy 



(216) 



s:^'=(c) = s'°'=(c) 



i 

1 



for arg(C) = , 



(217) 



s'°"(C) = s!!"(C) 



loc / 



for arg(C) = -7r/3 , 



(218) 



s^^(C) = s'°'=(c) 



for arg(C)=7r/3, and, 
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arg(Q=7i/3 



/ arg(Q=0 
arg(Q=-7i/3 



Figures. The contour T,ij^(^x-) in the -plane for Riemann-Hilbert Problem^. The boundary 
of the image ipiA{x){D), expanding as 0, is shown as a dashed curve. 



(219) 



for arg(C) = tt . 



4. Normalization: S'°'^(C) is normalized at infinity so that 



(220) 



s'°"(C)s'°"(C)"' 



as C — > cx) . 



We will describe how Riemann-Hilbert Problem ^ can be solved explicitly. First, we make an explicit 
change of variable to a new matrix T(C) by setting 



(221) 



si-(C) = T(C) 



g-(-C)-'/V2+W4 

g(-C)='/V2-»^/4 



If S'°'^(C) satisfies Riemann-Hilbert Problem ||, then it follows that T(C) obeys the following jump relations: 

-1 



(222) 



(223) 



(224) 



T+(C) = T_(C) 



1 



for arg(C) = , 



T+(C)-T_(C) 



1 1 
1 



for arg(^) = ±7r/3 , and, 



T+(C)=T_(C) 



1 
1 1 



for arg(C) = n . 



It is a general fact that NxN matrix functions that satisfy piecewise constant jump conditions, like the matrix 
T(C) does, can be expressed in terms of solutions of iVth order linear diiTerential equations with meromorphic 
(and often rational, or even polynomial) coefficients. In the 2x2 case, classical special functions therefore 
play a key role. In this case, we see immediately from the fact that the boundary values taken by S'°'^(C) on 



ASYMPTOTICS OF SEMICLASSICAL SOLITON ENSEMBLES: RIGOROUS JUSTIFICATION OF WKB 



37 



'^iA{x) s-re continuous that the matrix 



(225) Q(C):=— (C)T(C)-i 



JCloc o 

^(C)s'°^(C)-^-^(-C)^/^s'-(C) 



1 

-1 



gloc^^-)-l 



is analytic for C e C*. If we suppose that T(C) has a bounded derivative near the origin in each sector of 
C \ SiA(x) — a hypothesis that must be verified later — then we see that in fact Q(C) is analytic at the 
origin and is consequently an entire function of C. 

To work ou t ho w Q(C) behaves for large C, we need to use the normalization condition (22C) for S'°'^(C)- 
We interpret (|220| ) to mean both that 



(226) 

and also that 
(227) 



O 



s'°'=(c) 



loc 



-(0=1 + 



loc 



-(0 + 



s'°'=(C) . 



Both (22£) and (227) are again hyp othe ses that must be verified once we obtai n a s olution for S'°'^(C)- They 
are not true a priori by virtue of (220) alone; for example the decay rate in (22C) might not be as fast as 
and the error term might have rapid oscillations that would make its derivative larger than thus 
violating (227). It follows from our hypotheses that Q(C) must be a polynomial; in fact, 



(228) 



Q(C) 



1 



Under our hypotheses, it then follows that the matrix T(C) solves the linear differential equation 

-C 



(229) 



dC 4 



1 



T(C)- 



Upon introducing the new independent variable 



(230) 



2/3 



c 



we see that the elements of the second row of T satisfy Airy's equation: 



(231) 



and that the elements of the first row are given by 
(232) Tifc = - 



4\ dT2k 



d^ 



So which solutions of Airy's equation are the appropriate ones for our purposes? The first observation is 
that we need to specify different solutions of Airy's equation in each simply-connected region of the complex 
plane where 121(0 2^22(0 analytic. The assignment of solutions in these regions must be consistent 
with the jump conditions and asymptotics for T(0. From the jump conditions for T(0, we can see that in 
fact T2i(0 is analytic in C\M, while TjlLO is analytic except when arg(0 = ±27r/3 or ^ e M_. We now want 
to use the normalization condition ( P26[ ) to find sectors in which T21 and T22 are exponentially decaying. 
Then we will be able to uniquely identify these functions with particular solutions of Airy's equation that 
also decay. From the presumed asymptotic relation ( ^26[ ) we have (in terms of the variable Oi 



(233) 



7^21(0 



\ 1/6 2e 

32 ' 



V3 



1 + 



1 

0/2 
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and 

(234) 



1/6 



^1/4 



l + O 



1 



as ^ ^ oo with — tt < arg(^) < tt. These show that T2i(^) is exponentially decaying for — tt < arg(^) < — 7r/3 
and also for 7r/3 < arg(^) < tt, while T22(^) is exponentially decaying for —tt/3 < arg(^) < 7r/3. Significantly, 
both matrix elements are analytic throughout the sectors where they are exponentially decaying for large ^. 
As a basis of linearly independent solutions of Airy's equation we take the functions Ai(^) and Ai(^e^*'^/^). 
These decay in different sectors, and have the asymptotic expansions 



(235) Ai(0 = 
as ^ — > oo for — tt < arg(^) < tt, and 

(236) Ai(^e2*"/3) 



1 £-2? /3 

20F 

g-i7r/6 g2C3/2/3 



(1 + 0(1/0) 



-(1 + 0(1/0) 



as ^ ^ 00 for —n < arg(^) < —n/3. Comparing the expansion of Ai(^e^*'^/'^) with that of T2i(5) in the 
sector — TT < arg(^) < — 7r/3, we find that here 



(237) 



T2i(0 = -e-^"/i2gi/6^Ai(Ce2'-/3) . 



Since T21 is analytic in the lower half ^-plane, this relation holds identically for S>(^) < 0. Similarly, comparing 
the expansion of Ai(^) with that of T22(0 in the sector — 7r/3 < arg(^) < 7r/3, we find that here 



(238) 



T22(0 = e'^/%'/^V^Ai{0 . 



Being as T22 is analjd;ic for —2^/3 < arg(^) < 27r/3, this identity holds throughout the sector of analyticity. 
Restoring the original independent variable C, we find that for 9(^) > 0, 



(239) r2i(C) = -e-^/^^gi/e^^- 

and therefore throughout the same domain, 



2/3 



-i7r/3 



(240) 



1/6 



V^Ai' 



2/3 



t/3 



For C with 7r/3 < arg(C) < tt or — tt < arg(C) < — 7r/3, 

(241) T22(0 = e^'^/^e^/^'v^Ai 
and therefore throughout the same domain, 

(242) ri2(0=e^'"/'(f)'^'v^A: 



2/3 



2/3 



c 



The sector of C \ that is contained in both of these domains is 7r/3 < arg(^) < tt. It is sufficient 

to have specified the matrix elements of T{C) in this sector, since it may be consistently obtained in the 
remaining sectors of C \ by making use of the jump relations for T(^). The procedure is consistent 

because the cyclic product of these jump matrices is the identity: 



(243) 



" 1 


■ 




" 1 


1 ■ 


-1 


" -1 " 


-1 


" 1 


1 


1 


1 







1 




1 







1 



I. 
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Once T(C) is known for C G C \ the original unknown matrix S (C) solving Riemann-Hilbert 

Problem ^ is obtained directly by the transformation (221). It suffices to give a formula that holds for 
7r/3 < arg(C) < tt. We find 



(244) 



^lr(c) 

Sine) 

Sine) 



3 

32x1/6 
3 



1/6 



2/3 



— i7r/3 



2/3 



(-C)"'V2Ai 



2/3 



2/3 



c 



We will have found a solution to Riemann-Hilbert Problem ^ if we can verify the two hypotheses (c/. 
equations ( p26[ ) and ( p27| )) we made regarding the interpretation of the normalization condition (220) and 
the differentiability of S'°'^(C) at C = 0. One verifies these directly, using the explicit formulae given here. 

7.3.1. The local model for N(A) near A = iA{x). To build a better model for N(A) in the disk D than 
Nout(A), we begin by recalling the exact relationship between the matrix 0(A) and the matrix S(A): 



(245) 



0(A) = ^1^ S(^,^(,)(A))(7i 



for all X E D. The matrix N(A) is also explicitly related to 0(A). For A G I? in the part of the lens between 
the contours and Cm, 

Hn J 



1 —i exp 



(246) N(A) = 0(A)ctj; 

1 

for A G -D in the part of the lens between the contours Cm and 

'ie{\) 

J 11 exp I 

(247) N(A) = 0{X)a; ' 



tlN 







1 



and for all other A G I?, we simply have N(A) — 0(A). As we do not expect the difference between N(A) 
and the formal continuum limit approximation N(A) to be important in the disk D since it is isolated from 
the points A = and A = iA, we can obtain a guess for an approximation for N(A) that should be valid in 
D simply by substituting S{ipiA(x)W) for S{ifiA(x)(.^)) in these formulae. 

Putting these steps together, the model for N(A) for A g _D that we will use is defined as follows. For 
A e in the lens between and Ca/, set 



(248) Nendpoint(A) := (Ti = ^{^^A{x)W) 

for A G _D in the lens between Cm and L^;, set 

(249) Ne„dpoi„t(A) :=af^S((p,^(,)(A)) 



1 —I exp 


1 i exp 




Hn J 



1 
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and for all other X G D, set 



(250) Ncndpoint(A) CTi " S((^i^(^)(A))(Ti ■ 



1-J 



where S(C) is defined by ( pi5[ ). The most important properties of the matrix Nondpoint(A) in the disk D are 
the following. 

Proposition 16. The matrix Nendpoint(A) is piecewise analytic in the left and right half-disks of D. On the 
imaginary axis (which bisects D) oriented in the positive imaginary direction, 

(251) Ne„dpoint,-(A)-lNendpoi„t,+ (A) = N_(A)-1N+(A) , 



that is, the local model has exactly the same jump as N(A). For X Cz D, the matrix function Ncndpoint(A) is 



— 1/3 

bounded by a constant of order Hj^ . Also, for X £ dD, 



(252) Ne„dpoint(A)Nout(A)-i = I + 0(n]v^') . 



Proof: Computing the jump matrix for Ncndpoint(A) is straightforward. To show the bound on Nendpoint(A), 
we recall that ^'^°^{^PiA(x) {^)) is bounded by a quantity of order K^^^ , and note that the other factor of the 
same size comes from the factor S^°'^{LpiA(x){^)) via the normalization condition on this matrix and the fact 
that ^PiA{x){^) grows like hj^^^ for X £ D. Similar reasoning using (226) establishes the error in matching 
onto NoutC-^) on the boundary of D. □ 

8. The Parametrix and its Error 

We are now in a position to put all of our models together to build a guess for a uniformly valid approxi- 
mation of N(A). Such a guess is called a parametrix. 

8.1. Constructing the parametrix. To build the parametrix N(A) as a sectionally holomorphic matrix 
function, we simply combine the outer and local models. For all A satisfying |A| < where the parameter 
a is to be determined later, set 

(253) N(A) := NoHgin(A) . 
For A G _D, we set 

(254) N(A) := Nendpoi„t(A) , 
and by symmetry for all A <E D* we set 

(255) N(A) :=a2Nc„dpoint(A*)V2. 
Finally, for all remaining A e C, set 

(256) N(A) Nout(A) . 

The parametrix N(A) is holomorphic for A G C \ S, where S is the contour illustrated in Figure]^. 

8.2. Estimating the error. To determine the accuracy of the parametrix, we compare it directly with the 
original matrix N(A). That is, we consider the error matrix defined by 

(257) E(A) N(A)N(A)-i . 

This matrix is sectionally analytic in the complex A-plane, with discontinuities on a contour Y,e that is 
illustrated in Figure ^. Note that as a consequence of the symmetry of N(A) and N(A) under complex 
conjugation, we have E(A*) = (T2E(A)*(T2. If the parametrix is indeed a good model for N(A), then we must 
be able to show that the matrix E(A) is uniformly close to the identity matrix in the whole complex plane. 

While we do not know E(A) explicitly like we know N(A), we know from the normalization condition of 
both factors that 

(258) E(A) ^ I as X^ oo. 
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Figure 6. The contour Y,. The circles at the top and bottom of the figure are the boundaries 
of the disks D and D* respectively. The circle at the origin has radius h'^ . The contours 
to the left and right of the imaginary axis in the upper half-plane are portions of the lens 
boundaries Ll and Lr respectively. The remaining small segments present in the upper 
half-plane for |A| < ft^ are parts of Cl and Cr. 



It turns out that we can also calculate explicitly the ratio of boundary values taken by E(A) from both 
sides on each arc of Eg. That is, we know the jump matrix for E(A), and can express it explicitly in 
terms of N(A) and the jump matrix for N(A), both of which are knownQ. This means that the matrix E(A) 
itself is a solution of a particular Ricmann-Hilbert problem for which we know the data. By solving this 
Riemann-Hilbert problem, we will show that indeed E(A) is uniformly close to the identity matrix. 

There are two kinds of arcs in the contour E^: "matching" arcs of the circles dD, dD*, and |A| = h'^ 
where two different components of the parametrix have to match well onto each other, and the remaining 
arcs within the disks and outside the disks where the jump matrix for N(A) should be a good approximation 
to that of N(A). 

Consider one of the arcs of E^; oriented in some convenient way, and as usual let the subscript "+" 
(respectively "— ") denote a boundary value taken on the arc from its left (respectively right). We can easily 
see from the definition (b57|) that for A on this arc. 



(259) E+(A) =E_(A)ve(A) with ve(A) := N_(A)vn(A)vj;j(A)-1N_(A)-i 

where vn(A) and Vj;j(A) denote the jump matrices on the arc for N(A) and the parametrix N(A) respectively. 
If the arc under consideration is a "matching" arc, then the discontinuity in E(A) is wholly due to the 
mismatch of components of the parametrix, and the jump matrix vn(A) is therefore replaced with the 



identity matrix in (259). It then follows that an equivalent formula for ve(A) on a "matching" arc is the 
following: 

(260) ve(A) = N_(A)N+(A)"^ , for A on a "matching" arc of E^. 



^Or at least well- understood. We have characterized the parametrix for |A| < hf^ in terms of the matrix function -F(C) for 
which we have an existence proof and a characterization, but not an explicit formula. 
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Figure 7. The contour Tie- We have E C and the components o/ E^; \ E are shown in 
dashed lines to make a clear comparison with Figure |^. 

In this case, the two boundary values represent different components of the parametrix, for example Nout(A) 
would play the role of N+(A) and Nendpoint(A) would play that of N-(A) if the "matching" arc under 
consideration is an arc of dD, oriented clockwise. 
The key fact that we need now is the following. 

Proposition 17. The optimal value of the radius parameter a is a — 6/7. For this value of a, and for all 
v > arbitrarily small, 

(261) V£(A)-I = 0(4/'-'^). 
uniformly for all X £ Y,e- 

Proof: We begin by considering the "matching" arcs. We take the circle dD to be oriented in the clockwise 
direction. Here we find 

(262) ve(A) = Nendpoint(A)Nout(A)-i = I + 0{h]/') , 

with the error estimate coming from Proposition An estimate of the same form necessarily holds on the 
"matching" arcs of dD* according to the conjugation symmetry of E(A). The remaining "matching" arcs lie 
on the circle |A| = ft^, which again we take to be oriented in the clockwise direction. Here we find 

(263) ve(A) = Nongin(A)Nout(A)-i = 1 + 0(ft]f ") , 

with the error estimate coming from Proposition ^ 

We continue by considering the arcs of E^ with |A| < h'^. Using the fact recorded in Proposition Othat 



Norigin(A) has determinant one and is uniformly bounded, we see from (259) that the important quantity to 
estimate is simply vn(A)vj;j(A)^^ — I, the difference between the jump matrix ratio and the identity. First, 
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consider the portion of the contour Cl with |A| < 1. Using Proposition 13 and (|23|), we find that here 

1 



(264) vn(A)vj.(A) 



aL(A) 



r(l/2 + zC) 



1 



r(i/2-*c) 

where C = ipo(X). Using (156) and the uniform boundedness of the leading-order term on the right-hand 
side of (|156|) for |A| < /i^, we see that the matrix quotient in (264) differs from the identity matrix by an 



order ^ amount. VirtuaUy the same argument using (p5|) and (158) in conjunction with Proposition [iJ 



establishes that on the matrix quotient vn(A)vj;^(A) differs from the identity by a quantity of order 

Next, consider the contour Ll for |A| < hf^. On this contour, there is no jump for N(A), so 
I in the formula (|259| ) for ve(A). But from the formula for Vj;j(A) for A G Ll given in Proposition [T^, 



fc4a/3-l 

vn(A) 



we see by ordinary Taylor expansion that for |A| < hf^, we have ve(A) — I 



0(?i- 



""^) for A e ii. Virtually 



< K^. Using (|27| ) and the jump matrix Vp;f(A) for A € Cm recorded 



the same argument yields the same estimate for ve(A) — I on Ln with | A| < h%. Finally consider the contour 
Cm (the positive imaginary axis) with 
in Proposition |l^, we find that here 

(265) 



Vn(A)Vj;j(A) 



l + W+Z 



i0(X)/hMi 



W-)Z 



{W+ 

i{W+ + W-)aM{X)Z + W-Z 



e-'^W/hN (aM(A)-z(l + Z)) 
where Z is given in terms of ^ = Vo(A) by ( |18C| ) and where 
(266) W±:=l- e^^i+'^(0W-em/t^N _ 

Note that W+ and W- are both of order for |A| < hf^ by Taylor expansion arguments. Also, Z is 

uniformly bounded in the disk of radius K^, and e±'^(-^)/'»N both have modulus one for A G Cm in this disk. 
Also, from ( |16C1| ) we get aM{X) — «(! + Z) ^ 0{fi^^^ ^^'') for all > since 0(A) = on this part of Cj\/. 

ising from Taylor ; 

^.j-4a/3- 

o(V 



This error dominates those arising from Taylor approximation, and thus on Cm with |A| < h/^ we find that 

vn(A)vj^(A)~^ — I = 0(Sjlj^^^ ^ '^)- The corresponding estimates hold on the corresponding contours in the 
lower half-disk, by conjugation symmetry. Putting this information together with the uniform boundedness 
of Noiigin(A) and its inverse, we find that for all A G S^; with |A| < h'^, 



(267) 



ve(A) - I = 0(4"/'" 



for all 1/ > 0. 

Now we proceed to study the jump matrix ve(A) inside the disk D centered at the endpoint A = iA{x), 
assuming x ^ is fixed so that D is fixed and bounded away from the origin and from A = i A. The only 
contour we need to consider is the imaginary axis. The jump matrix ve(A) is given by ( |259| ). This time, 
the conjugating factors of Nondpoint,-(A) and its inverse are not uniformly bounded in D as Hn tends to 
zero. According to Proposition |l^ each conjugating matrix contributes an amplifying factor of hj^^'^ . Also 
according to Proposition |l^, we have that Vj<;j(A) is the same as the jump matrix for N(A). Therefore, using 
(H) and (H), we find that 



(268) 



Vn(A)Vj;j(A) 



(aA/(A) 



1 



Using Proposition ^ and the fact that 5R(0) < while e ^^W/^n jj^g^g modulus one on Cm within £), we get 
that VN(A)vj;f (A)~^ = I-\-0{h\^^) for all /it > 0. Combining this with the bounds on the conjugating factors, 
we find that within D, 



(269) 



ve(A)-I-0(4/^-'') 



for all /i > 0. By conjugation symmetry, the same estimate holds for the jump matrix ve(A) when A G D* . 
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Finally, we consider the parts of outside all of the disks, where we have set N(A) := Nout(A). On all 
of these parts of E^, Proposition ^ guarantees that the conjugating factors Nout,~(A) and Nout,-(A)~^ are 
uniformly bounded as tends to zero. So it remains to determine the magnitude of the difference between 
the quotient vn(A)vj;^(A)^^ and the identity. First consider the part of Cm between the disk at the origin 
and the disk D. Using Proposition ^to find Vj^(A) and recalling (p7|), we find 



(270) 



Vn(A)Vj;j(A) 



= while Q-^^W/fi-N jjas modulus one here then allows us to 
'') for all /X > 0. The a appears because we need the estimate 



Using Proposition ^ and the fact that ^(A) 
conclude that vn(A)V]<;j(A)^^ = I + 0{h]^°'' 
down to the outside boundary of the shrinking disk at the origin. Next we look at the contours Ll and Lr. 
On these contours there is no jump for N(A), and we see directly from Proposition |^ that vn(A)Vj;j(A)^^ — I 
is exponentially small as fi^r tends to zero through positive values. The next contour we examine is the 
portion of Cm lying above the disk D. Here we observe that Nout(A) has no jump, and because we have 
on this contour the strict inequality 3ff((^(A)) < 0, the matrix vn(A) is exponentially close to the identity 
matrix. To see this, note that here 



(271) 



vn(A) 



Am (A) 1 



because 6{X) = 0. While the relative error in replacing aM(A) by ie'^'^^^/^N jg j^q^ small near A = iA^ it is 
bounded. So the exponential decay afforded by the strict inequality on the real part of (/)(A) is maintained. 
Virtually the same arguments show that on the contours Cl and Cr outside of the disk at the origin, the 
quotient visi(A)vj<;j(A)~^ is again an exponentially small perturbation of the identity matrix, since on these 
contours there is again no jump of the parametrix Nout(A). Therefore, for all A e E^ outside all disks, we 
have 



(272) 



ve(A) 



0(4""-'^) 



for all /i > 0. 

We come up with an overall estimate for ve(A) 



for A e E^; by combining the estimates ( 262 ), ( 263 ) 



(267), ( p69D , and ( 272 ), and optimizing the error by choosing the parameter a. The optimal balance among 
all a e (3/4, 1) comes from taking a = 6/7, which gives an overall error estimate of 



(273) 



ve(A) - 



uniformly for all A G E^;, for all v > Q. This proves the proposition. □ 

The following is then a consequence of the L/^ theory of Riemann-Hilbert problems (see the analogous 
discussion in pMMOl). 



Proposition 18. For Km sufficiently small, the Riemann-Hilbert problem for E(A) has a unique solution. 
Let R > be sufficiently large so that T,e is contained in the circle of radius R centered at the origin. Then 
uniformly for all A outside of this circle, and for all ^ > however small, the matrix E(A) satisfies 

(274) E(A) - I = Oih]!''") 

with the .size of the matrix measured in any matrix norm. 

Proof: Only one thing must be verified in order to deduce existence and uniqueness from the general 
theory: the Cauchy-kernel singular integral operators defined on the contour E^;, which depends on Hn 
because of the shrinking boundary of the circle at the origin, have L'^(T,e) norms that can be bounded 
uniformly in N. But this fact follows in this case from the fact that for N sufficiently large, the circle 
|A| = h'^ intersects only radial straight-line segments (we chose the contours to all be exactly straight lines 
in some fixed neighborhood of the origin), so that the portion of E^; near the origin simply scales with h"^. 
The uniform bound we need can then be established using the fact that the Cauchy operators commute with 
scaling. A similar result was established under more general conditions in |KMMOC|. Once existence and 
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uniqueness have been estabhshed, the estimate of E(A) follows from an integral representation formula for 
this matrix: 

(275) E(A)=I+-^/ (s-A)-i(m(s)(vE(s)-I))ds 

•J Y^E 

in which the matrix function m(A) for A G is an element of L^(I]£;) with a norm that is bounded 
independently of hj^ (and in fact converges to the identity matrix in l^iJ^E) hj^ tends to zero). □ 

We are now in a position to prove our main result, which we presented as Theorem |l| in the introduction. 

Proof of Theorem^: We have been setting i = all along, so the function ip^" (x) given by ( |l5|) can be found 
from the matrix N(A) by the relation (cf. equation (^)) 

(276) ^o"ix) = 2i lim AiVi2(A). 
Writing N(A) = E(A)N(A) N(A) + (E(A) - I)N(A), we get 

(277) ^o^ix) = 2i lim AiVi2(A) + 2i lim A(£;il(A) - l)7Vi2(A) + 2i lim Ai;i2(A)iV22(A) . 

A — *oo A — >oo A — >oo 



Now, the first term on the right-hand side of ( p77D can be evaluated explicitly since near A = cxd we have 
N(A) = Nout(A) = 0(A), and we have an explicit formula {cf. equation (^) for 0(A). We find 

(278) 2i lim AiVi2(A) = A{x) , 

A — 'oo 

which is the "true" initial data that we started with, before making any modifications based on the WKB 



approximation of the spectral data. When we consider the second term on the right-hand side of ( P77D , we 
see that as a consequence of the normalization of the matrices E(A) and N(A), 

(279) E(A) =I-hO(l/A), and N(A) = I + 0(1/A) , 



as A — > oo. Therefore the second term on the right-hand side of (E77^ vanishes identically. Finally, for the 



third term on the right-hand side of (277) we can again apply the normalization condition for N(A) to obtain 

(280) 2i lim Ai;i2(A)iV22(A) = 2i lim A£;i2(A). 

A — s-oo A — >oc 

Putting these steps together, we have 

(281) i^o^ix) - A{x) = 2i lim A£;i2(A) . 

A — >oo 

The proof of the theorem is finished upon using the integral formula ( ^75| ) for E(A) and Proposition |lj. □ 

9. Discussion 

Using the new technique of simultaneous interpolation of residues by two different analytic interpolat- 
ing functions, combined with "steepest-descent" techniques for matrix-valued Riemann-Hilbert problems, 
we have established the validity of the formal WKB approximation of the spectrum in the nonself adjoint 
Zakharov-Shabat eigenvalue problem (|^) in the sense of pointwise convergence of the potentials. Strictly 
speaking, our analysis applies to certain classes of potential functions whose most important property for 
our purposes is their real analyticity, and then we obtain convergence for all nonzero values of x. 

In order to extend the result of Theorem |l| to a; = 0, some different steps are required. Since A{x) — s- A 
as a; ^ 0, the local analysis that we carried out independently for A ~ (cf. § |7.2|) and for A ~ iA{x) {cf. 



7.3) will need to be combined. Consequently, a different local model for N(A) will need to be constructed 



near X = iA = iA{0). Due to the presence of the gamma functions in the asymptotics established in § 7.2 
and given in Proposition |l^ and Proposition |l^, it is likely that the construction of the local model will 
require knowledge of the solution of a new Riemann-Hilbert problem that, like that for the matrix F{C) in 
§7.1, cannot be solved explicitly. Nonetheless, one expects that to establish the validity of Theorem |l] for 
X = will require only technical modifications of what we have done here. 

Understanding the nature of the WKB approximation at the level of the potentials is one step in a larger 
ongoing program to obtain corresponding information at the level of the (unknown) spectrum itself. Indeed, 
quantifying the difference between the true spectrum of a given potential A{x) and the WKB approximation 
of the spectrum, in terms of motion of eigenvalues, will be necessary before it can be proven that the rigorous 



46 



P. D. MILLER 



asymptotic analysis of SSEs is relevant to the problem of semiclassical asymptotics for the initial-value 
problem for the focusing nonlinear Schrodinger equation (|l|) . One imagines that a study of the semiclassical 
limit for should proceed by first generating from the given initial data ip{x, 0) — A{x) the corresponding 
well-defined SSE, and then using the fact that by combining the results of | KMMO(| with Theorem]^ from 



this paper, one has a complete picture of the limiting behavior of the SSE for an open interval of time t that 
is independent of H and includes t = 0, and moreover that according to Theorem ^ the SSE is pointwise 
close to the given initial data A{x) for t — 0. The problem here is that the focusing nonlinear Schrodinger 
equation is known to have modulational instabilities whose exponential growth rates become arbitrarily large 
in the semiclassical limit. There is, therefore, the very real possibility that while the SSE is close to the 
initial data A{x) at t = 0, it is not close to the corresponding solution of (|l|) for any positive t. In order to 
control the difference for positive time, it is necessary to know in advance how much the SSE spectral data 
differs from the true (unknown) spectral data, as it is the spectral data that is the starting point for analysis 
( cf. Riemann-Hilbert Problem |l|) . One can imagine obtaining this sort of information from the pointwise 
estimate given in Theorem |] by Rayleigh-Schrodinger perturbation theory. 
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